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Die Erfindung betrifft ein Verfahren zum Bestimmen 
der Position eines Erdpunktes nach dem Oberbegriff 
des Patentanspruches. Es wird mit Hilfe von Funksigna- 5 
len, die von erdumkreisenden Satelliten gesendet war- 
den, durchgefuhrt Derartige Satelliten senden unter- 
schiedlichste Signale. Am einfachsten auszuwerten sind 
Signale, bei denen der Tr^ger mitgesendet wird und bei 
denen der Modulationskode bekannt ist. Am schwierig- 10 
sten ist dagegen die Standortbesiimmung unter Ausnut- 
zung von Funksignalen, bei denen der Trager unier- 
driickt ist und bei denen die Modulationsart unbekannt 
ist. Derartige Signale sind 2. B. die» die im NAVSTAR 
Global Positioning System (GPS) gesendet werden. Bei 15 
diesen Signalen sLndert sich der Modulationskode in 
nicht bekannter Weise. 

Ein Stand der Technik» von dem das eingangs genann- 
te Verfahren ausgeht, ist in Bull. Geod. 53, 1979, Seiten 
1 1 7 — 1 38, beschrieben. 20 

Einige Systeme zur Standortbestimmung durch Funk 
machen Gebrauch von der Richtcharakteristik des 
Strahlungsdiagramms einer Sende* oder Empfangsan- 
tenne. Andere Systeme, einschlie&lich die Erfindung, 
stQtzen sich nicht auf die Richtcharakteristik irgendei- 25 
ner Antenne. Die Erfindung geh5rt zu der allgemeinen 
ICategorie von Systemen, bei denen der Standort einer 
Empfangsantenne durch Messen des Unterschieds zwi- 
schen den Phasen oder den Gruppenlaufzeiten oder bei- 
den von Signalen bestimmt werden kann. die von zwei 30 
oder mehr unterschiedlichen Sendeantennen. dessen 
Standort schon bekannt ist, eintreffen. Falls zwei Sende- 
quellen synch ronisiert sind oder falls die Synchronisa- 
tionsabweichung der zwei Sender unabh^ngig bekannt 
ist, wird dann durch Messung des Gruppenlaufzeitun- 15 
terschieds der von den beiden Quellen ankommenden 
Signale der Empfanger tokalisiert, und zwar dreidimen- 
sional auf einem besonderen Rotations-Hyperboloid, 
dessen Brennpunkte die Standorte der Sender darstel- 
len. Falls ahnliche, am gleicheri Empfangsort vorgenom- 40 
mene Messungen von Signalen mehrerer unterschiedli- 
cher, geeignet aufgestellter Sender kombiniert werden, 
kann der Empfangsort eindeutig anhand des Schnitt- 
punktes der entsprechenden Hyperbolioden bestimmt 
werden. 45 

Techniken zur Bestimmung der relativen Lage von 
unterschiedlichen Orten, und zwar einer in bezug auf 
den anderen, durch Messung der Phasen- oder der 
Gruppenlaufzeitunterschiede zwischen Funksignalen, 
die gleichzeitig an diesen Orten empfangen werden, sind 50 
ebenso im Stand der Technik bekannt und werden ins- 
gesamt als Techniken der Geodasie durch Funkinterfe- 
rometrie erwahnt. Die Antennen an den einzelnen Or- 
ten bilden ein Interferometer, und der relative Lagevek- 
tor, der sich von einer Antenne zu der anderen erstreckt, 55 
wird als Basisvektor des Interferometers bezeichnet. 
Der Basis- oder Relativlagevektor zwischen zwei An- 
tennen kann gewohnlich mit geringerer Unsicherheit als 
die Lage jeder einzelnen Antenne bestimmt werden, da 
viele potentielie Fehlerquellen die Neigung zeigen, auf eo 
die Messungen an beiden Antennen fast gleich einzu- 
wirken und demzufolge, wenn die Differenz zwischen 
beiden Antennen gebildei wird, sich aufheben. Die Geo- 
dasietechnik mittels Mikrowellenfunkinterferometrie 
ist bekannt fiir ihre unerreichte Kombination aus Ge- 65 
nauigkeii. Geschwindigkeit und Reichweite zur Bestim- 
mung von Relativlage- oder Interferometerbasisvekto- 
ren. Eine derartige Bestimmung kann entweder auf 



Messungen des Grundla^^Bunterschieds oder auf 
Messungen des PhasenunteWffieds, oder beiden Unter- 
schieden der Signale gegrundet werden. die an den bei- 
den Enden des Basisvektors empfangen werden. Pha- 
senmessungen sind von Natur aus genauer als Gruppen- 
laufzeitmessungen, jedoch ist die Interpretation der 
Phasenmessungen weitaus schwieriger infolge ihrer in- 
newohnenden, ganzperiodischen Mehrdeutigkelt. Eine 
allgemelne Dtskussion (iber InterferometermeQtechni- 
ken und den damit verbundenen Problemen der Inter- 
pretation wird in einem Artikel "Radio Astrometry^, in 
Annual Reviews of Astronomy and Astrophysics, Vol. 
14 (1976), pp. 197—214, von Charles C. Counselman HI 
diskutiert. Eine grofle Sammlung von relavanten techni- 
schen Abhandlungen erscheint in Conference Publica- 
tion 2115 der National Aeronautics und Space Admini- 
stration, mit dem Titel "Radio Interferometry Techni- 
ques for Geodesy^. Geodasie durch Radiointerferome- 
trie wurde mit Radiosignalen praktiziert, die von ver- 
schiedenen Quellen einschlieBlich naturlicher, wie z. B. 
Quasare, und kunstlicher, wie z. B, Satelliten des NAV- 
STAR Global Positioning System (GPS), stammen. 

Wie bekannt, gibt es gegenwartig etwa sechs erdum- 
kreisende GPS-Satelliten. Die Umlaufbahnen der Satel- 
liten konnen mit einer Genaugkeit von etwa 2 Metern 
bestimmt werden. Diese Satelliten emittieren Funksi- 
gnale mit Wellenlsingen in der Nslhe von 19,0 Zentime- 
tem und ebenso 24,4 Zentimetem. Geht man davon aus, 
daB die ganzperiodischen Mehrdeutigkeiten der interfe- 
rometrischen Phasenbeobachtungen dieser Signale ein- 
wandfrei ausgelost werden, kann der Basisvektor, der 
sich von einer Antenne zur anderen erstreckt, interfero- 
metrisch mit einem Unsicherheitsfaktor bestimmt wer- 
den, der viel kleiner als die Wellenlangen der GPS- 
Obertragungen ist Durch Bestimmung von drei Basisli- 
nien, wobei jede Basis eine LSnge in der GroQenord- 
nung von 100 Metern aufweist, mit Hilfe der interfero- 
metrischen Phasenmessungen von GPS-Sigrialen wurde 
gezeigt, daB man eine Genauigkeit innerhalb etwa 1 
Zentimeters erzielen kann^ und zwar gemaQ eines Be- 
richtes, veroffentlicht in Eos (Transaction of the Ameri- 
can Geophysical Union), Vol. 62, Seite 260, 28. April 
1981, von Charles C. Counselman HI, S. A. Gourevitch, 
R. W. King, T. A. Herring, 1. 1. Shapiro, R. L. Greenspan. 
A. E. E. Rogers, A. R. Whitney, und R. J. Cappallo. Das 
bei diesen interferpmetrischen Basisbestimmungen ver- 
wendeten Verfahren griindete auf der bekannten Tech- 
nik der direkten, an einer zentralen Stelle vorgenomme- 
nen Kreuzkorrelation der Signale, die getrennt, aber 
gleichzeitig an den beiden Enden jeder Basis empfangen 
werden. 

In der US-PS 41 70 776 wird ein System zum Messen 
von Anderungen eines Basisvektors zwischen zwei Or- 
ten auf der Erde unter Verwendung von Signalen be- 
schrieben. die von den GPS-Satelliten gesendet werden, 
bei dem die an jedem Ort empfangenen Funksignale 
genau zeitmarkiert und anschlieBend uber Teiefonlei- 
tungen zu einer zentralen Stelle ubertragen werden, wo 
ein naher Echtzeit-Phasenvergleich durch Kreuzkorre- 
lation beider Signalreihen durchgefuhrt wird. Das in die- 
sem Patent beschriebene System beinhaltet Parabolre- 
flektor-Empfangsantennen. Da die FunkfluBdichte eines 
GPS-Signals relativ zu dem HintergrundrauscHpegel 
gering ist und da die Bandbreite eines GPS-Signals bei 
weitem die Bandbreite einer Telefonleitung uberschrei- 
tet, ist das Signal-Rauschverhaltnis der uber die Tele- 
fonleitung von jedem Ort ubertragene Leistung gering. 
Es ist daher zur Anhebung dieses Signal-Rauschverhalt- 
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nisses auf einen sinnvollen P^^wesentlich, daB Para- 
bolantennen mit groQen SamWmachen bei diesem Sy- 
stem verwendet werden. Ein anderer wichtiger Grund 
fur die Verwendung derartiger Antennen ist darin zu 
sehen, dafi diese eine Richtwirkung aufweisen» so daB 5 
Signale, die die Antennen anders als direkt von der gcr 
wunschten Quelle erreichen, abgewiesen werden. 

Systeme zum Messen von Basisvektoren. die andere 
Arlen von Signalen von erdumkreisenden Satelliten 
verwenden, sind ebenso bekannt. 10 

In einem Artikel ''Miniature Interferometer Termi- 
nals for Earth Surveying** (MITES), in Bulletin Geodesi- 
que. Vol. 53 (1979), Seiten 139-163, von Charles C 
Counselman III und Irwin I. Shapiro wird ein vorge- 
schlagenes System zum Messen von Basisvektoren be- 15 
schrieben, das vielfrequente Funksignale verwendet, die 
von erdumkreisenden Satelliten gesendet werden wtir- 
den, wobei in diesem System die Phasen der empfange- 
nen Signale getrennt an jedem Ende der Basis bestimmt 
werden. D. h. das an einer Stelle empfangene Signal 20 
wird nicht kreuzkorreliert mit dem an der anderen Stel- 
le empfangenen Signal, um die Phasendifferenz zwi- 
schen den beiden Signalen zjii bestimmen. Um die Pha- 
senmehrdeutigkett zu t5sen, stQtzt sich das MITES-Sy- 
stem auf die Kombination von Messungen an einen Satz 25 
von bis zu zehn Frequenzen, die geeignet zwischen 1 
und 2 GHz verteilt sind. Jedoch gibt es, soweit bekannt, 
keine gegenwartig die Erde umkreisende Satelliten, die 
derartige Signale aussenden. 

Systeme zum Messen der Relativlage unter Verwen- 30 
dung von Signalen, die von anderen Quellen als kunstli- 
chen Satelliten gesendet werden,;sind ebenso bekannt. 
Ein Beispiel fOr ein derartiges System, das eine auf den 
Mond basierende Obertragung verwendet, ist ebenso in 
der US-PS 41 70 776 offenbart. 35 

Systeme zum Messen entweder einer einzelnen Lage 
Oder einer Relativlage unter Verwendung von Signalen 
von Quellen, die nicht von erdumkreisenden Satelliten 
stammen, sind ebenso bekannt. Z. B. wird in einem Arti- 
kel von W, O. Henry, "Some Developments in Loran" in 40 
journal of Geophysical Research, Vol. 65. Seiten 
506—513, Febr. 1960, ein System zur Bestimmung einer 
Lage bzw. eines Standorts (z. B. eines Schiffes auf See) 
beschrieben, und zwar unter Verwendung von Signalen, 
die von auf dem Erdboden stationierten Sendern stam- 45 
men. Dieses System, bekannt als Loran-C- Navigations- 
system, weist mehrere tausend Kilometer lange Ketten 
von synchronisierten Sendern auf, die auf der Erdober- 
flache stationiert sind, wobei alle Sender die gleiche 
Tragerfrequenz, 100 Kilohertz, verwenden, und wobei 50 
jeder Sender durch ein eindeutiges, periodisches Im- 
pulsmuster amplitudenmoduliert wird. Dieses Muster, 
das Vorzeichenumkehrungen der Amplitude ein- 
schlieQt, ermoglicht dem Empfanger, die Signale von 
verschiedenen Sendern zu unterscheiden, Eine geeigne- 55 
te Kombination von Beobachtungen von mehr als ei- 
nem Paar von Sendern kann eine Bestimmung der Emp- 
fangerlage auf der Erdoberflache lief ern. 

Ein anderes Beispiel eines derartigen Systems stellt 
das Omegasystem dar, das in einem Artikel yon Pierce eo 
mit dem Titel *'Omega", in IEEE Transactions on Aero- 
space and Electronic Systems, Vol. AES-l, no. 3, Seiten 
206— 215, Dez. 1965, beschrieben wird. Beim Omegasy- 
stem werden die Phasendifferenzen der empfangenen 
Signale im Prinzip fast wie die Gruppenlaufzbiten. in 65 
dem Loran-OSystem gemessen. Da die verwendeten 
Frequenzen sowohl in dem Loran-C- und dem Omega- 
system sehr hiedrig sind, sind die Genauigkieiten der 
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Lagemessungen mit died^^Vstemen im Vergleich zu 
den erwahnten Satellitens^TOmen ziemlich schlecht. 

Der Stand der Technik beinhaltet auch andere Ver- 
fahren zur Bestimmung des Standorts und der relativen 
Lage mit Hilfe des Globalposiiionssystems (GPS). Das 
Standardverfahren, beschrieben z. B. in einem Artikel 
von J. J. Spitker, Jr, in einem Artikel in Navigation, Vol, 
25, no. 2 (1978), Seiten 121 - 146. und weiterhin beschrie- 
ben in mehreren anderen in der gleichen Ausgabe dieses 
Journals erscheinenden Ariikeln, basiert auf Messungen 
der Unterschiede zwischen den Gruppenlaufzeiten oder 
den "Zeiten" des Empfangs der codierten Modulation 
der GPS-Signale. Im Prinzip stellt dieses Verfahren ein 
hyperbolisches Positionsverfahren dar und ist im we- 
sentlichen ahnlich dem von Loran. Die Bandbreite von 
etwa 10 MHz der GPS- Modulation begrenzt die Ge- 
nauigkeit der Gruppenlaufzeitmessun^ und folglich der 
Lagebestimmung durch das Standardverfahren auf 
mehrere Dezimeter. Eine Genauigkeit in der GroBen- 
ordnung von einem Zentimeter ist moglicherweise 
durch Verwendung von Tragerphasenmessungen er- 
reichbar, wie dies z. B. in einem Artikel von ]. D. Bossier, 
C. M. Goad und P. L Bender mit dem Titel "Using the 
Global Positioning System for Geodetic Positioning", in 
Bulletin Geodesique, Vol 54. no. 4, Seite 553 (1980). 
beschrieben wird Jedoch hat jedes veroffentliehte Ver- 
fahren den Nachteil. daB die Kenntnis und die Verwen- 
dung der Code modulation, die chiffriert werden kann, 
erforderlich ist, oder daB eine Kreuzkorrelation der an 
verschiedenen Orten empfangenen Signale erforderlich 
ist. oder daB die Verwendung groBer Antennen erfor- 
derlich ist, um den empfangenen Rausch- bzw. Storab- 
stand anzuheben und die Interferenz von reflektierten 
Signalen zu unterdrticken, oder das Verfahren weist 
mehr als einen dieser Nachteile auf. 

Der Erfmdung liegt die Aufgabe zugrunde, ein Ver- 
fahren zur Standortbestimmung mittels GPS-Signalen 
anzugeben, also Funksignalen, bei denen der Tr^ger un^ 
terdruckt und der Modulationskode nicht bekannt ist. 

Das erfindungsgemaBe Verfahren ist diirch die Merk- 
male des Patentanspruchs gegeben. 

Von besonderer Bedeutung fur das erfindungsgema- 
Be Verfahren ist, daB es mit einer einzigen Antenne 
gleichzeitig die Signale von mehreren Satelliten mi&t, 
daB es fQr diese Signale Vorhersagen fur Doppelver- 
schiebungen macht und daB es die Signale mit sich selbst 
korreliert, was vorzugsweise durch Korrelieren der Sei- 
tenbander eines jeweiligen Signals erfolgt. 

Das erfindungsgemaBe Verfahren weist keinen der 
Nachteile bekannter Verfahren zum Auswerten von 
GPS-Signalen auf. d. h. es muB keine groBe oder stark 
bundelnde Empfangsantenne verwendet werden, das 
am^Standort empfangene Signal muB nicht mit einem an 
einer anderen Stelle empfangenen Signal kreuzkorre- 
liert werden, und es ist keine Kenntnis des Kodes erfor- 
derlich, die die Tragersignale modulieren. 

Die Erfmdung wird im folgendeh anhand der Zeich- 
nungen naher erlautert. Es zeigt: 

Fig. 1 eine Vorrichtung zur Bestirnmung eines Basis- 
vektors durch Funkinterferometrie mittels GPS-Satelli- 
ten gemaB der Erfindung; 

Fig. 2 ein Blockdiagramm eines der in Fig. 1 gezeig- 
ten Interferometer- Feldterminals; 

Fig. 3 ein Blockdiagramm der in Fig. 2 gezeigten An- 
tennenanordnung; 

Fig. 4 ein Blockdiagramm der in Fig. 2 gezeigten 
Empfangereinheit; 
Fig. 5 ein Blockdiagramm der in Fig. 2 gezeigten Di- 
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gitalelektronikeinheit; 

Fig. 6 ein Btockdiagramm de^PPig. 5 gezeigten Si- 
gnalutnsetzers; 

Fig. 7 ein Blockdiagramm eines Korrelatormoduls 
der in Fig. 5 gezeigten Korrelationsanordnung; 5 

Fig. 8 ein Blockdiagramm eines numerischen Oszilla- 
tormoduls der in Fig. 5 gezeigten numerischen Oszilla- 
toranordnung und 

Fig. 9 ein Blockdiagramm des in Fig. 2 gezeigten 
Feldterminal-Rechners. 10 

Die Erfindung ist auf eine Technik zum Messen des 
Basisvektors zwischen zwei Punkten. wie z. B. Vermes- 
sungsmarken, auf der Erde durch Funkinterferometrie 
unter Venvendung von Doppelseitenbandfunksignalen 
mit unterdrOcktem Trager, die von erdumkreisenden 15 
Satelliten des NAVSTAR Globaipositiorissystems 
(GPS) stammen, gerichtet. Diese Technik schlieBt das 
Messen der Phasen der Tragerwellen, die in den an je- 
dem Ort empfangenen Signalen enthatten sind, und an- 
schlieBend das Bearbeiten der Phaseninformation, die 20 
an beiden Stellen erhalten wird, zur Bestimmung des 
Basisvektors ein. Ein Vorteil dieser Technik besteht dar- 
in, daS man die Tr^gerphasen ohne Bezug auf die 
Kenntnis der codierten Signale miBt, die in den Satelli- 
ten zur Modulation der Trager verwendet werden. Ein 25 
anderer Vorteil besteht darin» da& es nicht erforderlich 
ist, die empfangenen Signale, entweder in Echtzeit oder 
durch Transport der Aufzeichnungen» von den zwei Or- 
ten zu einem gemeinsamen Ort zu transportieren. Ein 
weiterer Vorteil besteht darin, daB die Verwendung von 30 
groBen oder eine starke Richtwirkung aufweisenden 
Antennen nicht erforderlich ist Noch ein weiterer Vor- 
teil besteht darin, daB diese Technik relativ immun ge- 
gen Fehler ist, die durch Streuung und Reflexionen der 
Funkwellen in der Nahe der Empfangsantennen hervor- 35 
gerufen werden. 

Obwohl die Erfmdung anschlieBend spezieli ftir die 
Verwendung mit GPS-Satelliten beschrieben wird, ist 
verstandiich, daB gewisse Aspekte hiervon nicht allein 
auf die Verwendung mit derartigen Satelliten begrenzt 40 
sind und in Verbindung mit Signalen, die von anderen 
Quellen empfangen werden, gut verwendbar sein kdn- 
nen. 

Wie bekannt, umkreisen Satelliten des NAVSTAR 
Globalpositionssystems (GPS) die Erde in einer Hohe 45 
von etwa 20 000 Kilometern und ubertragen Signale in 
einem Frequenzband, das bei 1575,42 MHz zentriert ist, 
bekannt als das "LI "-Band, sowie Signale in einem zwei- 
ten Band, das bei 1227,60 MHz zentriert ist, bekannt als 
das "1-2"- Band. Die Signale werden derart moduliert, 50 
daB nahezu symmetrische Ober- und Unterseitenbander 
erzeugt werden, wobei der Trager vollstandig unter- 
druckt wird. 

FOr jedes Band kann das Signal von einem bestimm- 
ten Satelliten, das an einer bestimmten Stelle empfan- 55 
gen wird. als eine Funktion der Zeit in folgender Form 
dargestellt werden: 



s(t) = m(t)cos(2 Ttfot -f O) + n(t) sin (2 nfot -h <D) 
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wobei m(t) und n(t) Modulationsfunktionen darstellen, 
und zwar jeweils^eine reelwertige Zeitfunktion: fo stellt 
die Tragersollfrequenz dar, die gleich 1575,42 MHz fur 
LI und 1227,60 MHz fur das L2-Band ist; (I) stellt die 
empfangene Tragerphase in Radian dar, die unbekannt 65 
ist und bestimmt werden soli. Jede der Modulations- 
funktionen m(t) und n(t) stellt eine pseudofallige Zeit- 
funktion mit einem Miitelwert von Null dar; Die beiden 



W 

Funktionen sind gegenseiti^^Btwinklig bzw. orthogo- 
nal, jede der Funktionen^K zur Modulation des 
Ll-Tragers fur irgendeinen Satelliten verwendet wird, 
ist auch rechtwinkiig zu der entsprechenden Funktion. 
die fur jeden anderen Satelliten verwendet wird, obwohl 
fur einen bestimmten Satelliten die gleiche m|[t)- oder 
n(t)-Funktion oder beide verwendet werden konnen, um 
sowohl die LI - und die L2-Trager zu modulieren. 

Die Bandbreiten der beiden Funktionen m(t) und n(t) 
unterscheiden sich durch eihen Faktor von genau 10, 
wobei m(t) die engere und n(t) die breitere Bandbreiie 
aufweist. Gewohnlich sind bei LI sowohl m(t)- und 
n(t)-Signalkomponenten vorhanden, und bei L2 ist nur 
die n(t)-Komponente vorhanden, die m(t)-Funktion wird 
auf Null gesetzt oder "ausgeschaltet**. Die Leistungs- 
spektraldichte von m(t), die dem Modulationssignal ent- 
spricht, das in der GPS-Literatur als "clear/acquisi- 
tion^-Code bekannt ist, ist proportional zu der Funktion 

sin^(7tF/l,023MHz) 
(kF/1,023 MHzy 



wobei F die Modulationsfreiquenz darstellt Diese Funk- 
tion hat eine Halbwertsbreite bei halbem Maximum von 
etwa 450 kHz. D. h. der Funktionswert ist etwa 0,5 fur 
F=- ±450 kHz, wohingegen der Wert fur F=0 Eins ist 
Die Leistungsspektraldichte von n(tX die dem Modula- 
tionssignal entspricht, das in der GPS-Literatur als "pre- 
cise code** oder "p code** bekannt ist, ist proportional zu 

sin^ (7iF/10,23MHz) 
(71F/IO.23 MHz)2 



Somit betr^gt die Halbwertbreite bei halbem Maxi- 
mum der Leistungsspektraldichte von n(t) etwa 
4,5 MHz. 

Fiir das LI -Signal von 1575,42 MHz ist der mittlere 
Quadratwert von n(t) gewohnlich gleich ein Halb dem 
von m(t);d. h. 

<n^{t)> ^ 0,5 <m2(t)>. - 

(Es ist mdglich, einen GPS-Satelliten in auBerordentli- 
chen Betriebsarten zu betreiben, bei denen das Verhaltr 
nis der mittleren Quadratwerte oder das Leistungsver- 
haltnis sich von Wert 0,5 unterscheidet. insbesondere ist 
ein Wert von Null moglich,) Demzufolge entspricht das 
Verhaltnis der Leistungsspektraldichte von n(t) zu der 
von m(t) gewohnlich etwa 0,5 10 = 0,05 fiir einen Wert 
von F nahe Null, so daB, falls ein an das Spektrum m(t) 
angepaBtes BandpaBfilter bei der Ll-Tragerfrequenz 
zentriert ist, etwa 90 Prozent der Leistung, die in dem 
Ausgang dieses Filters enthalten ist, von der m(t)-Si- 
gnalkomponente stammen wird und weniger als 10 Pro- 
zent von der n(t)-Komponente. Zur Vereinfachung wird 
in der ubrigen Beschreibung davon ausgegangen. daB 
das GPS-Ll -Signal keine n(t)-Komponente aufweist 
und die folgende einfachere Form aufweist: 

s(t) « m(t) cos (2 Tifot -I- <D). 

Im allgemeinen ist, die empfangene Tragerphase <t> 
eine sich langsam andernde Zeitfunktion, so daB die 
tatsachlich empfangene Tragerfrequenz durch die fol- 
gende algebraische Summe gegeben ist: 
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f = fo + (2n)-^(d<D/dt), 



wobei fo die Tragersollfrequenz und dO/dt die Zeitablei- 
tung von 0 darstellt. Durch den Ausdruck ''sich langsam 
ahdernde" ist gemeinti daQ (2n)"* (dO/di) sehr klein im 5 
Vergleich zu fo und zur Bandbreile von m(t) ist. Der 
Hauptgrund ftir die zeitliche Anderung von d> ist in der 
Dopplerverschiebung zu sehen, wodurch f von fo bis zu 
etwa plus Oder minus 4»5 kHz abweichen kann. 

Das empfangene Signal s(t) enthalt keine diskrete to 
spektrale Leistungskomponente bei der Tragerfre- 
quenz, da der Mittelwert von m(t) Null ist. Soniit wird 
der Trager vollstandig unterdriickt, und die Leistungs- 
spektraldichtefunktion des LI -Signals s(t) ist gleich der 
Leistungsspektraldichtefunktion der Modulation m(t), 15 
ubertragen vom Basisband auf die empfangene Trager- 
frequenz f. Da m(t) eine reellwertige Zeitfunktion dar- 
stellt, ist dessen Leistungsspektraldichte eine gerade 
symmetrische Funktioh der Frequenz. Somit hat die Lei- 
stungsspektraldichte von s(t) eine gerade Symmetric be- 20 
ziiglich der Tragerfrequenz f und wird als Doppelseiten- 
bandspektrum bezeichnet. Der Teil dieses Leistungs- 
spektrums, der den Frequenzen grofier als f entspricht, 
wird oberes Seitenband genannt; der Teil, der den unte- 
ren Frequenzen entspricht, wird als un teres Seitenband 25 
bezeichnet. (Die leichte Aymmetrie, bestenfalls etwa 3 
Telle von 10^ zwischen dem oberen und dem unteren 
Seitenband infolge der Doppler-**Dehnung*' des Signals 
ist hier nicht von Bedeutung.) 

GemaQ der Erfindung wird eine Antenne an jedem 30 
Ende eines Basisvektors aufgestellt. Die Signale, die 
durch jede Antenne empfangen werden» werden in obe- 
re und untere Seitenbandkbmponenten getrennt. Diese 
getrennten Komponenten werden gefiltert, in Ein-Bit- 
Digitalform umgewandelt und anschHeBend zusammen 35 
multipliziert Ihr Produkt wird digital durch Korrelation 
mit Quadraturausgangen eines Emfanger- bzw. Oberla- 
gerungsoszillators analysiert, urn die Leistung und die 
zu diesem Uberlagerungsoszillator relative Phase der 
Tragerwelle zu bestimmen, die in dem von jedem Satel- 40 
liten empfangenen Doppelseitenbandsignal enthalten 
ist Unterschiede in der Dopplerverschiebung. werden 
verwendet, um die Trager der unterschiedlichen Satelli- 
ten zu unterscheiden. Somit werden die Leistungen und 
Tragerphasen der Signale mehrerer Satelliten gleichzei- 45 
tig gemessen, und an jeder Vermessungsmarke werden 
numerische Daten erhalten, die die MeBergebnisse re- 
prasentieren. Die Messungen werden bei jeder Marke 
in Echtzeit durchgefiihrt, und zwar phne Bezug auf die 
Signale, die an irgendeinem anderen Ort empfangen 50 
werden und ohne Kenntnis irgendeines der codierten 
Signale. die die GPS-Trager modulieren. Die Daten der 
gleichzeitig, aber unabhangig an den beiden Vermes- 
sungsmarken durchgefuhrten Messungen — und zwar 
einmal pro Sekunde fur eine Zeitspanne von ausreichen- 55 
der Dauer. wie z, B, von etwa 5000 Sekunden — werden 
dann zusammen verarbeitet, um den sich von einer Mar- 
ke zu anderen erstreckenden Basisvektor zu bestimmen. 
Zwei Verarbeitungsverfahren werden offenbart Bei 
dem eirien Verfahren wird eine "Mehrdcutigkeitsfunk- eo 
tion" berechnet, die eine Funktion der MeBdaten und 
eines Ausgangs- bzw. Versuphsvektors B des Basisvek- 
tors ist. Der Vektorraum von 6 wird systematrsch abge- 
sucht, um den eindeutigen Wert von & zu finden, der die 
berechnete Funktion maximiert. Dieser Wert von 6 wird 65 
als erwunschte Bestimmung des unbekahnten Basisvek- 
tors bgenommen. 

Mit Bezug auf Fig. 1 wird nun eine Vorrichlung 11 . 
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zum Bestimmen eines B^^Hctors b^emaQ der Erfin- 
dung verdeutlicht. Der ^^Kvektor b, der manchmal 
auch anschlieOend als "Basislinie" bezeichnet wird, stellt 
den Relativlagevektor von einer, Vermessungsmarke 
SM-1 in bezug auf die andere Vermessungsmarke SM-2 
dar. Die Basis erstreckt sich von der Vermessungsmarke 
SM-1, die sich an dem Ursprung oder einem Ende der 
Basislinie befindet, zu der Vermessungsmarke SM-2, die 
sich an dem Endpunki oder dem anderen Ende der Ba- 
sislinie befindet Die Vorrichtung It umfaBtzwei rntelli- 
gente Interferometer-Feldterminals 13-1 und 13-2, und 
zwar eines an jedem Ende der Basislinie. sowie einen 
Rechner, der strukturell und funktionell in eines der 
Terminals 13 eingebaut und Teil davon sein kann oder 
der, wie gezeigt eine getrennte Einheit 15 darstellen 
kann. 

Die Vorrichtung erfordert zum normalen Betrieb ge- 
wisse numerische Daten von externen Quellen und 
ebenso einige Einrichtungen zur Obertragung numeri- 
scher Daten zwischen dem Rechner 15 und jedem Ter- 
minal 13, und zwar vor und nach oder (fakultativ) wah- 
rend der Durchfuhrung der Basismessungen. 

Ehe mit den Messungen zur Bestimmung der Basis; 
begonnen wird, werden Daten von einem ersten Daten- 
speicher 17, und zwar steJlvertretend fur die Umlauf- 
bahnen mehrerer GPS-SatelHten — von denen zwei mit 
GPS-1 und GPS-2 bezeichnet zu Illustraiionszwecken 
gezeigt werden — in den Rechner 15 eingegeben, und 
zwar zusammen mit Approximationsdaten, die stellver- 
tretend fur die Orte der Vermessungsmarken SM- 1 und 
SM-2 sind. und die von einem zweiten Datenspeicher 19 
erhalten werden. Die letzteren Daten konnen z. B. die 
Vermessungsmarkenorte innerhalb einer Genauigkeit 
von einigen Kilometern reprasentieren. Von diesen Sa- 
tellitenbahn- und MeQortdaten erzeugt der Rechner 15, 
in Tabellenform als Funktion der Zeit, eine Voraussage 
der Dopplerfrequenzverschiebung, die das von jedem 
GPS-Satelliten gesendete 1575.42-MHz-Signal beim 
Empfang an jeder Vermessungsmarke haben wird. Der 
Rechner 15 erzeugt ebenso eine tabellenformige Vor- 
aussage des Leistungspegels des Signals, das von jedem 
Satelliten an jeder Marke empfangen wird. Die voraus- 
gesagte Leistung ist Null, falls der Satellit sich unterhalb 
des Horizonts befindet Die vorausgesagte Leistung 
stellt eine Funktion des vorausgesagten Hohenwinkels 
des Satelliten uber dem Horizont dar infolge der Win- 
kelabhangigkeit des Gewinns einer Empfangsantenne 
(bei der Vermessungsmarke) und, gewdhnlich in gerin- 
gerem AusmaO. einer Sendeantenne (am Satelliten). Die 
Tabellen mit den vorausgesagten Frequenzverschie- 
bungen und -leistungen werden — fur eine Zeitspanne, 
die die Zeit fiir die erwarteien Messungen umfal^t — fur 
alJe GPS-Saielliten, die nach aller Erwartung bei jeder 
MeOmarke sichtbar sind, durch irgendeine bekannte 
Einrichtung, wie z. B. durch eine Telefon- oder Funkte- 
lefonzwischenleitung zu einem Speicher eines kleineren 
Rechnerij ubertragen und eingegeben, der in dem be- 
sonderen Interferometer-Feldterminal 13, das an dieser 
Vermessungsmarke aufgestellt wird oder bereits aufge- 
stellt wurde, enthalten ist; Alternativ konnen die Fre- 
quenz- und Leistungsvorhersagetabellen durch den im 
Interferometer-Terminal enthaltenen. Rechner erzeugt 
werden. 

Die Dopplerfrequenzvorhersagen werden nach For- 
meln berechnet, die im Stand der Technik gut bekannt 
sind. Die FehlergroBen bei solchen Vorhersagen sind in 
der GroQenordnurig von 1 Hertz pro Fehlerkilometer in 
bezug aiif den angenommenen Ort der Vermessungs- 
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marke. Der zus^tzliche Fehler ^^^r Frequenzvorher- 
sage infolge eliies Fehlers bei^^ Extrapolation der 
SateDitenbahn liegt normalerweise in der GroBenord- 
nung von 1 Hertz oder weniger fiir wenigstens einen 
Tag vorher gemachte Vorhersagen. Frequenzvorhersa- 5 
gefehler bis zu mehreren Hertz sind im Zusammenhang 
mit der Erfindung zul^ssig. Die Vorhersagen uber die 
empfangene Leistung brauchen nicht sehr genau zu 
sein; Fehler von mehreren Dezibel sind zulilssig* da die- 
se Voraussagen nicht fiir einen sehr kritischen Zweck 10 
verwendet werden. Sie dienen in der Hauptsache dem 
Feldterminalrechner dazu, zu priifen, ob das gewunsch- 
te Signal, und nicht irgendein Storsignal empfangen 
wird. Unter mdglicherweise einem gewissen Verzicht an 
Zuverlassigkeit konnten die Leistungsvorhersagetabel- 15 
len eliminiert werden. 

Ein InteHferometer-Feldterminal 13, das an einer Ver- 
messungsmarke aufgestellt wurde. empfangt nun gleich* 
zeitig die 1575,42-MHz-SignaIe von einer Vielzahl von 
Satelliten, und zwar bis zu sieben, aber keinesfalls weni* 20 
ger als zwei Satelliten. Fiir eine genaue Bestimmung der 
zu erlangenden Basis ist es fiir die Terminals an beiden 
Enden der Basis wesentiich, die Satelliten gleichlaufend 
zu beobachten. 

Elektronikschaltkreise (nachfolgend beschrieben) in- 25 
nerhalb jedes Terminals trennen die empfangenen Si- 
gnaie in obere und untere Seitenbandkomponenten auf 
und analysieren unter Verwendung der Vorhersage der 
Dopplerfrequenzverschiebung diese Seitenbandkom- 
ponenten, um die Leistung und die Phase der Trager- 30 
welle zu bestimmen, die in dem von jedem Satelliten 
empfangenen Signal enthalten ist. Die Daten von diesen 
Leistungs- und Phasenbestimmungen werden innerhalb 
des Feldterminals gespeichert und evtl. durch irgendei- 
ne bekannte Einrichtung an den zentralen Rechner 15 35 
wieder zurQckgegeben. 

Die Daten von den beiden Interferometer-Feldtermi- 
nals 13-1 und 13-2 mussen zusammen verarbeitet wer- 
den, um eine genaue Bestimmung des Basisvektors zu 
erhalten, 40 

Es ist zu bemerken, daB fiir den Betrieb dieses Sy- 
stems keine Einrichtung fur eine Ferniibermittlung oder 
Obertragung der Daten erforderlich ist Die Terminals 
13-1 und 13-2 konnen kdrperlich zu der gleichen Stelle 
wie der Rechner 15 transportiert werden, und dort kon- 45 
nen die vorhergesagten Tabellen vom Rechner 15 zu 
den Terminals 13 Qbertragen werden. AnschlieBend 
konnen die Terminals 13, die die Tabellen in ihren Spei- 
chern enthalten, zu den Vermessungsmarken SM-1 und 
SM-2 befdrdert werden, wo die Satelliten beobachtet 50 
werden. Nach Beendigung dieser Beobachtungen kon- 
nen die Terminals 13 zu der Stelle des Rechners 15 
zurQckbefdrdert werden, wo die Tr^gerphasendaten 
von beiden Terminals zu dem Rechner zur Bearbeitung 
ubertragen werden kdnnen, 55 

Mit Bezug auf Fig. 2 werden die Hauptkomponenten 
eines Interferometerterminals 13, auch als "Feldtermi- 
nar bezeichnet. verdeutlicht. Jedes Feldterminal 13 
weist eine Antennenanordnung 21 auf, die mit einer 
Elektronikanordnung 23 mit Hilfe eines Koaxialkabels 60 
25 verbunden ist. 

Jede Antennenanordnung 21 weist eine Antenne 27 
und eine Vorverstarkeranordnung 29 auf. Die Antenne 
ist auf der Vermessungsmarke SM positioniert und kann 
in, Art der Antenne, die in der US-Patentanmeldung 55 
3 23 328, eingereicht am 20. November 1981, beschrie- 
ben ist, aufgebaut sein. Ohne Rucksicht auf die Art der 
IConstruktion muQ der Ort des Phasehzentrums die An- 
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tenne 27 beziiglich der Ve^^Bungsmarke SM genau 
bekannt sein. Die Antenne^^Tin dieser Patentanmel- 
dung beschrieben ist, ist in dieser Hinsicht ausreichend 
Die Ungenauigkeit in der Positionierung ihres Phasen- 
zentrums betragt bestenfalts einige Millimeter. 

Die Antenne 27 empfSngt die 1575.42- MHz- Funksi- 
gnale, die von den GPS-Sateliiten gesendet werden. Die 
empfangenen Signale werden durch den Vorverstarker 
29 verstarkt und iiber das Koaxialkabel einer Empfan- 
gereinheit 31 zugefuhrt, die in der Elektronikanordnung 
23 enthalten ist, wobei die Empfangereinheit 31 einen 
Seitenbandseparator 33, einen Empfangerarbeitsstrom- 
kreis 34 und einen Oszillatorschaltkreis 35 einscl^liefit. 

In dem Seitenbandseparator 33 wird der obere Sei- 
tenbandteil der Signale, die denjenigen Teil der von al- 
ien kombinierten Satelliten empfangenen Signale ent- 
halten. der einen Bereich von Funkfrequenzen innehat, 
der sich von 1575.42 MHz nach oben erstreckt, von dem 
unteren Seitenbandteil getrennt, der Funkfrequenzen 
unterhalb 1575,42 MHz entspricht. Zur Durchfuhrung 
dieser Trennung verwendet der Seitenbandseparator 33 
ein 1575,42-MHz-BezugsssignaI, das von dem Oszilla- 
torschaltkreis 35 geliefert wird. 

Die Empfangereinheit 31 liefert drei Signale in ainalo- 
ger Form zu einer Digitalelektronikeinheit 37. Ein ana- 
loges Signal, bezeichnet mit u(t), reprasentiert die obere 
Seitenbandkomponente der empfangenen Funkfre- 
quenzsignale, die ins Basisband iibertragen wurde. Das 
zweite analoge Signal, gekennzeichnet mit e(t), repra- 
sentiert die untere Seitenbandkomponente, die ebenso 
ins Basisband iibertragen wurde. ]edes dieser beiden 
Signale enthalt Beitrage von alien sichtbaren Satelliten. 
Das dritte, an die Digitalelektronikeinheit 37 gelieferte 
Signal ist ein sinusfdrmiges Signal mit einer Frequenz 
von 5,115 MHz» das die AusgangsgroBe eines freisch- 
wingenden, stabilen Quarzkristalloszillators in dem Os- 
zillatorschaltkreis 35 darstellt. Die AusgangsgroBe die- 
ses gleichen Oszillator wird hinsichtlich der Frequenz 
mit einem festen ganzen Faktor von 308 innerhalb der 
Oszillatoranordnung multipliziert, um die Bezugsfre- 
quenz von 1575,42 MHz zu erhalten, die von dem Sei- 
tenbandseparator verwendet wii-d. Die Genauigkeit der 
Frequenzen, die durch die Oszillatoranordnung 35 er- 
zeugt werden, betragt in der Regel etwa ein Teil von 10^ 
obwohl eine Genauigkeit von ein Teil von 10^ zul^ssig 
ware. 

In der Digitalelektronikeinheit 37 wird jede der drei 
analogen EingangsgroBen in ein digital-logisches Signal 
ubergefiihrt Die Digitalsignale werden unter der Steue- 
rung eines Feldterminalrechners 39 verarbeitet, um die 
Tragerleistungs- und Phasendaten zu erzeugen. Die Di- 
gitalelektronikeinheit 37 ist mit dem Feldterminalrech- 
ner 39 mit Hilfe eines bidirektionalen Datenbusses 41 
verbunden. Der Feldterminalrechner 39 kann ein Digital 
Equipment Corporation (DEC)-Modell LSI— 11/2 Mi- 
crocomputer sein; in diesem Fall kann der Datenbus 41 
ein DEC-"Q"-Bus sein. 

Die Tragerphasendaten werden in dem Speicher des 
Feldterminalrechners 39 gespeichert, bis diese Daten zu 
dem zentralen Rechner 15 zur Verarbeitung iibertragen 
werden. Wie bereits bemerkt, kann der zentrale Rech- 
ner 15 eliminiert werden und die Verarbeitung in einem 
der Feldterminalrechner 39 durchgefuhrt werden. Die 
Phasendaten konnen ebenso durch den Feldterminal- 
rechner 39 auf ein Datenspeichermedium, wie z. B. eine 
Magnetbandkassette oder eine Magnetplatte (nicht dar- 
gestellt), herausgeschrieben werden. Die Daten konnen 
ebenso uber eine direkte elektrische Verbindung oder 
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Ilfonverbindung oder 
bhtungen ubertragen 



uber ein Modem und einc 
durch viele andere ubliche 
werden: 

Mit Bezug auf Fig. 3 werden weitere Einzelheiten der 
Komponenten der Antennenanordnung 21 gezeigt. Die 5 
Antennenanordnung 21 beinhattet eine Antenne 27, die 
— wie erwahnt — derart aufgebaut ist» daB ihr Phasen- 
zentrum genau beziiglich der Vernnessungsmarke posi- 
tioniert werden kann. Die 1575;42-MHz-Funksignaie, 
die durch die Antenne 27 empfangen werdea werden 10 
der Vorverstarkerschaltung 29 zugefiihrt, dessen Funk- 
tion es ist, ihren Leistungspegel ausreichend anzuheben, 
um die Dampfung des Koaxialkabels 25, das die Anten- 
nenanordnung 21 mit der Empfangereinheit 31 verbin- 
det, und das Hintergrundrauschen, das in dem Eingangs- 15 
verstSrker in der Empfangereinheit 31 erzeugt wird, zu 
beseitigen. 

In der Vorverstarkerschaltung 29 werden die von der 
Antenne 27 empfangenen Signale zuerst durch ein 
BandpaOf liter 43 von etwa 50 MHz Bandbreite, zen- 20 
triert bei 1575,42 MHz, gefiltert Die Funktion des Fil- 
ters 43 besteht darin, eine Oberlastung der Empfanger- 
einheit 31 durch starke Storsignale zu verhindern, die 
auBerhalb des GPS-Signalbandes vorhanden sein kon- 
nen. Die AusgangsgroBe des BandpaBfilters 43 wird ei- 25 
nem passiven Diodenbegrenzer 45 zugefOhrt, der zum 
Schutze des Durchbrennens eines rauscharmen Verstar- 
kers 47 dient, und zwar infolge von irgendweichen, sehr 
starken Signaien, wie z. B. Signale, die von in der Nahe 
befindlichen Hochleistungsradaranlagen abgestrahlt 30 
werden konnten, Der rauscharme Verstarker 47 stellt 
einen normalen Gallium-arsenid-Feldeffekttransistor 
(FET)- Verstarker mit einer Rauschzahl von etwa 2 db 
dar. 

Die Gleichstromleistung fur den rauscharmen Ver- 35 
starker wird tiber das Koaxialkabel 25, das mit der Vor- 
verstarkerschaltung 29 verbunden ist, von der Empfan- 
gereinheit 31 zugefuhrt, und zwar uber eine Funkfre- 
quenzdrossel 49 und einen Spannungsregler 51. Ein 
Kondensator 53 koppelt die Funkfrequenzausgangsgro- 40 
Be des rauscharmen Verstarkers 47 mit dem ICabel 25, 
wahrend es den Gleitstrom von dem Verstarker abhalt. 

Mit fiezug auf Fig. 4 werden detaillierter die Kompo- 
nenten der Empfangereinheit 31 dargestellt. Die Emp- 
fangereinheit 31 schlieBt einen Empfangerarbeitsstrom- 45 
kreis 34, einen Seitenbandseparator 33 und einen Oszil- 
latorschaltkreis 35 ein. Der Empfangerarbeitsstrom- 
kreis liefert Gleichstrom fur den Betrieb des Osziliator- 
schaltkreises 35, des Seitenbandseparators 33 und uber 
das Koaxialkabel 25 zum Betrieb des rauscharmen Ver- 50 
starkers 47 in der Antennenanordnung 21. Der Oszllla- 
torschaltkreis 35 liefert eine Bezugsfrequenz von 
1575.42 MHz zu dem Seitenbandseparator 33 sowie ei- 
ne Bezugsfrequenz von 5,115 MHz aii die Digitalelek- 
tronikeinheit 37. Der Seitenbandseparator 33 trennt die 55 
Signale, die in einem Funkfrequenzband empfangen 
werden, das bei 1575,42 MHz zentriert ist und sich von 
dieser Frequertz nach oben und unten erstreckt, in ge- 
trennte obere und iintere ^ieitenbandkomponenten im 
Basisband - 60 

Der Empfangerarbeitsstromkreis34 enthalt geregelte 
Gleichspannungsnetzgerate 61 und zusatzlich eine Ak- 
kumulatorbatterie63. Die Batterie 63 ermdglicht es, das 
dem Kristalloszillator 65 in dem Oszillatorschaitkreis 35, 
der Echtzeituhr in der Digitalelektronikeinheit 37 und 65 
dem Datenspeicher des Feldterminalreehners 39 ohne 
Unterbrechung Leistung zugefuhrt werden kann, und 
zwar ungeachtet von Unterbrechungeii des externen 
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Hauptstromversorgungs^^B, die auftreten kdnnen; 
somit wird die Frequenzsraffitat des Oszillators erhal- 
ten, und die Uhrzeiteinstellung sowie die in dem Rech- 
nerspeicher gespeicherien Daten gehen nicht verloren. 

Der Oszillaior 65 in dem Oszillatorschaitkreis 35 stellt 
einen Quarzkristallbszillator dar, wie z. B, ein Frequen- 
cy and Time Systems (FTS)-Modell 1001 dar, das eine 
Ausgangsfrequenz von 5.115 MHz innerhalb eines 
Bruchteils von 10^ oder weniger liefert. Das FTS-Modell 
1001 besitzt eine Stabilitai von etwa einem Teil von 10*** 
pro Tag und einem Teil von 1 0* ^ uber Zeitintervalle von 
1 bis 100 Sekunden, und ist demzufolge fur diese An- 
wendung mehr als geeignet 

Der Oszillator 65 liefert zwei identische Ausgangs- 
groBen, und zwar eine an die Digitalelektronikeinheit 37 
und die andere an einen t575,42-MHz-Normalfrequenz- 
generator 67 in dem OszillatorschaUkreis 35. 

Der 1 575,42- MHz-Normalfrequenzgenerator 67 ent- 
halt einen spannungsgesteuerten Transistoroszillator 
(VCO) 69, der bei einer Frequenz von 393,855 MHz 
schwingt; d. h. 77 x 5,1 15 MHz. Die Phase des Oszilla- 
tors wird bezuglich der Phase des 5,1 15-MHz-Bezugssi- 
gnals durch die Wirkung einer Phasenverkettungs- 
schleife stabilisiert, die aus dem VCO 69, einem Koppler 
71, einer Divisionsschaltung 73, einem Phasenfrequenz- 
fehlerdetektor 75 und einem Schleifenfilter 77 besteht. 
Ein Teil der Ausgangsleistung des VCO 69 wird mil 
Hilfe des Kopplers 71 an deh Eingang der Frequenzdivi- 
sionsschaltung bzw. Frequenzteilerschaltung 73 gekop- 
pelt,die aus standardisierten emittergekoppelten Logik- 
IC-Schaltungen besteht. die eine Division durch 11 und 
anschlieBend durch 7 vornehmen. Der Ausgang der Fre- 
quenzteilerschaltung 73 stellt ^die "variable** Eingangs- 
gr5Be und der 5,1 15- MHz- Ausgang des Oszillators 65 
stellt den "Bezugs^-Emgang fur den standardisierten 
ECL-IC-Pha$enfrequenzdetektor 75, wie z. B. Motoro- 
la-Typennummer MC 12 040, dar. Der Ausgang des De- 
tektors 75 wird in dem Schleifenfilter 77 tiefpaBgefiltert, 
um die Steuerspannung zu erhalten, die die Eingangs- 
groBe des VCO 69 darstellt Der Ausgang des VCO 69 
wird vervierfacht hinsichtlich der Frequenz durch eine 
Folge von zwei standardisierten, symmetrischen Dio- 
denfrequenzverdopplern 79 und durch einen Verstarker 
81 verstarkt, um die 1575,42-MHz-Ausgangsfrequenz zu 
erhalten, die den Seitenbandseparator 33 ansteuert 

Die in einem bei 1575,42 MHz zentrierten Frequenz- 
band enthaltenen Signale, die uber das Koaxialkabel 25 
von der Antennenanordnung 21 an dem Eingang 83 des 
Seitenbandseparators 33 empfangen werden, werden 
durch einen Gleichstrom-Sperrkondensator 85 an ein 
BandpaBfilter 87 angekoppeit und mitlels eines Ein- 
gangsverstarkers 89 verstarkt. Die Gleichstromleistung 
fur den VoryerstSrker 29 (in der Antennenanordnung) 
wird iiber eine Funkfrequenzdrossel 91 und dem Emp- 
fangerarbeitsstromkreis 34 auf das Koaxialkabel 25 ge- 
koppelt 

Der HF-Leistungsverteiler oder die"Verzweigung" 
93, die 1575,42-MHz-Oberlagerungsoszillator-Quadrat- 
urgabelschaltung 95, die beiden Ringmischer 97 und 99 
und die Breitbandvideofrequenzquadraturgabelschal- 
tung 101 in dem Seitenbandseparator 33 umfassen einen 
dualen Funkfrequenz-Basisbandumformer fur ein Ein- 
seitenband oder einen ^^Demodulator" vom bekannten 
"Phasenabgleich"-Typ. Ein derartiger Demodulator 
wurde beispielsweise in einem Artikel von Alan E. E. 
Rogers in Proceedings of the IEEE, vol. 59 (1971),Seiten 
1617 — 1618 beschrieben. Seine Arbeitsweise wird hier 
nachfolgend beschrieben. 
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Man bezeichnet mit fo die ^^Aenz des Bezugssi- 
gnals, das.durch den Oszillator^ffaltkreis zu dem Sei- 
tenbandseparator 33 geliefert wird. Die Frequenz fo ent- 
spricht nominell 1575,42 MHz, was der Tragersollfre- 
quenz der GPS-Satelliten-"Lr-Obertragungen ent- 5 
spricht, und zwar vor der Dopplerverschiebung (erster 
Ordnung). Oann konnen die AusgangsgroBen 102 und 
103 der Quadraturgabelschaltung 95 als sin 2 nfot bzw. 
cos 2 Tifpt geschrieben werden. Die AusgangsgroBen, die 
um SO*" phasenverschoben sind, stellen die "Oberlage- 10 
rungsoszillator (local oscillator)**- EingangsgrdBen fur 
die Mischer 97 bzw. 99 dar. Die Hochfrequenzeingangs- 
groBen fiir die beiden Mischer sind identisch. Die Basis- 
bandausgangsgroBen der Mischer sind demzufolge 
identisch mit Ausnahme einer Phasenverschiebung von 15 
n/2 Radian. (Mit "Basisband'* ist der naher bei Null, al$ 
bei fo gelegene Frequenzbereich gemeint, der dem Un- 
terschied zwischen der Eingangsfrequenz und fo ent- 
spricht) Die Richtung dieser Phasenverschiebung, Pha- 
senvoreilung oder -nacheilung, ist davon abhangig, ob 20 
die Eingangssignalfrequenz oberhalb oder unterhalb 
von fo liegt Somit ist es moglich, entweder die Obersei- 
tenband-(£igenfrequenz hdher) oder die Unterseiten- 
bandeingdnge auszuwSlhlen und das gegentiberliegende 
Seitenband durch Verschiebung der Phase eines Mi- 25 
scherausgangs um zusatzlich n/2 Radian auszuscheiden 
und dann die beiden Mischerausgange entweder zu ad- 
dieren oder zu subtrahieren (abhangig davon, welches 
Seitenband gewunscht ist). 

Die Quadratur- bzw. 90''-Phasenverschiebungsgabel- 30 
schaltung 101, die zwei Eing^nge 109 und lit sowie 
zwei Ausgange 105 und 107 hat, fiihrt diese n/2-Phasen- 
verschiebung und diese Addition/Subtraktion durch. 
Der obere Ausgang 105 der Gabeischaltung 101 ist ge: 
geben durch die arithmetische Summe des oberen Ein- 35 
gangs 109 plus des unteren Eingangs 111, wobei beide 
Eingange in der Phase um einen von der Frequenz ab- 
hangigen Betrag verzogert wurden, jedoch mit einer in 
bezug auf den oberen Eingang, um konstant nJ2 Radian 
groBeren Phasenverschiebung des unteren Eingangs. 40 
und zwar abhangig von der Frequenz. Der untere Aus- 
gang 107 ist durch die arithmetische Differenz der glei- 
chen beiden, unterschiedlich phasenverschobenen Ein- 
gange 109 und 111 gegeben, wobei die Differenz im 
Sinne von oben minus unten genommen wird. Der fest- 45 
gesetzte, itll Radian (eine Viertel Periode), Phasenun- 
terschied wird genau fiir alle Frequenzen zwischen fnp 
und wenigstens fLP eingehalten, wobei fnp— 10 kHz viel 
geringer als ftp ^450 kHz und ftp etwa gleich der ein- 
seitigen Bandbreite der GPS-"C/A "-Modulation m(t) ist, 50 
wie vorstehend diskutiert wurde. Der Aufbau einer 
Quadraturgabelschaltung, die diese Eigenschaften auf- 
weist, ist in dem zitierten Artikel von Rogers dargestelk. 

Nun werden die Ausgange der Quadraturgabelschal- 
tung 101 getrennt durch identische Videoverstarker 113 55 
und 1 15 verstarkt, durch HochpaBfilter 117 und 119 so- 
wie TiefpaBfilter 121 und 123 gefiltert Die Filter 117 
und 1 19 sind identische HochpaBfilter mit einer Nieder- 
frequenzsperre bei (hp- DerZweck dieser HochpaBfilter 
117 und 119 besteht darin, Gleichstromkomponenten 60 
und irgendwelche Niederfrequenzspektralkomponen- 
ten der Mischerausgange zu eliminieren, und zwar be- 
zuglich Frequenzen, die annahernd gleich oder niedri- 
ger als die maximal mogliche Grofle der Dopplerver- 
schiebung, die ein GPS-Sateilitsignal haben kann, sind. 65 

Es ist erwunscht. derartige Komponenten auszu- 
schlieBen, da andererseits diese die nachfolgende Be- 
stimmung der empfangenen, dopplerverschobenen Tra- ' 
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gerphase in der Digitalelek^^Keinheit und dem Rech- 
ner des Feldterminals behln^ffn konnten. Solche poten- 
tiellen Interferenz- bzw. Storsignale k5nnen niederfre- 
quentes Funkelrauschen enthalten, das in den Mischern 
selbst erzeugt wird, oder konnen von einer ICombina- 
tion aus einer Unausgeglichenheit der Mischer und (un- 
erwunschten) niederfrequenten Amplituden- oder Pha- 
senschwankungen des 1 575.42-MHz-Bezugssignals 
oder des Verstarkungsfaktors irgendwelcher Funkfre- 
quenzsignalverstarker herruhren, die den Mischern vor- 
angehen. Eine andere potentielle Quelle einer niederfre- 
quenten Interfrequenz bzw. Storung stellt das Brummen 
Oder die Welligkeit der Ausgangsspannungen bzw. 
-strome des Netzgerates dar. Eine weitere Quelle konn- 
te ein storendes Dauer (CW)-Signal sein. dessen Fre- 
quenz Nahe bei fo Hegt. 

Die TiefpaBfilter 121 und 123 stellen identische Tief- 
paBfilter mit einer Bandbreite gleich fLP, gleich der ein- 
seitigen Bandbreite von m(t), dar. Das Ansprechverhal- 
ten jedes Filters, als Funktion der Frequenz, ist derart 
zugeschnitten, daB diese mit der Leistungsspektraldich- 
te von m(t) zusammenpaBt Der Zweck dieser Filter 
besteht darin, Rauschen und Interferenz auBerhalb der 
Bandbreite von m(t) auszuschlieBen. Es ist zu bemerken, 
daB GPS-T code^-Modulationssignale n(t) mit groBer 
Bandbreite hier normalerweise eine Interferenz- bzw. 
St5rquelle darstellen wurden. Zumeist wird etwa 80 
Prozent der von dem Signal n(t) stammenden Leistung 
durch diese TiefpaBfilter abgehalten. Dieser Abwei- 
sungsgrad ist ausreichend, um sicherzustellen, daB die "P 
code'*-Interferenz eine vernachlassigbare Wirkung hat. 
Wir mochten jedoch bemerken, daB fur den Fall. daB die 
SchmaIband-m(t)-Modulation in den GPS-Satelliten ab- 
geschaltet wurde, dann die Breitband-n(t)-Modulation 
nicht linger ein unerwiinschtes Stdrsignal darstellen 
wiirde; es wurde das gewiinschte Signal werden. Ein 
derartiger Schalter konnte in der GPS-Signalstruktur 
durch VergroBerung der Bandbreiten der TiefpaBfilter 
35 um einen Faktor von 10 untergebracht werden, um 
diese an das neue "Signal** anzupassen. 

Die AusgangsgroBe u(t)des TiefpaBfilters 121 repra- 
sentiert die nach unten umgesetzte und gefilterte obere 
Seitenbandkomponente des Ursprungssignals s(t). Die 
AusgangsgroBe l(t) des TiefpaBfilters 123 reprasentiert 
das untere Seitenband. Es sollte bemerkt werden, daB 
das Spektrum von u(t) in bezug auf die Frequenz nach 
oben verschoben wird und das Spektrum von l(t) bezug- 
lich der Frequenz nach unten verschoben wird, und 
zwar relativ zu dem Spektrum der ursprunglichen Mo- 
dulation m(t) um einen Betrag gleich (f— fo), d. h. dem 
Unterschied zwischen der aktuell empfangenen Trager- 
frequenz f und der Oberlagerungsoszitlatorfrequenz fo. 
[Falls die Dopplerverschiebung des Tragers (f— fo) ne- 
gativ ist, wird dann das u(t)-Spektrum nach unten und 
das l(t)-Spektrum nach oben verschoben.] Man nimmi 
an, daB die GroBe dieser Verschiebung kleiner als fHP 
und sehr viel kleiner als ftp sein soil. Dieser Annahme 
wird genugt, falls die Frequenzverschiebung primar von 
der Dopplerverschiebung herriihrt, die niemals die Gro- 
Be von 5 Kilohertz Obersphreiten kann, vorausgesetzt, 
daB fHP etwa gleich 10 kHz gesetzt wird. Irgendeine 
Absetzung der Frequenz des Bezugskristalloszillators 
65 von der gewiinschten 5,1 15-MHz-Frequenz wird 
auch eine (308mal groBere) Verschiebung der u(t)- und 
l(t)-Spektra hervorrufen. Normalerweise wird jedoch 
eine derartige Verschiebung sehr viel kleiner als fHP 
sein. 

Zusatzlich zu der Frequenzverschiebung der oberen 
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und unteren Seitenbandausg^^pu(t) und l(t) entsteht 
eine frequenzabhangige, dis^^ierende Phasenver- 
schiebung jedes Ausgangs infolge der Quadraturgabel- 
schaltung 101. Jedoch ist fiir den speziellen Quadrature 
gabelschaltungsaufbau von Rogers (oben zitiert) diese 5 
Phasenverschiebung zu gefing, um wichdg zu sein, Ahn- 
lich warden die zus^tzlichen Phasenverschiebungen, die 
durch das BandpaBfilter 87 und die Hoch- und TtefpaB- 
filter 117, 119, 121 und 123 hervorgerufen werden, unbe- 
deutend sein. falls Standardfiiterkonstruktionen ver- 10 
wendet werden. Jeder dieser Effekte neigt zur.Aufhe- 
bung, wenn die Differenz zwischen den Terminals in der 
nachfolgenden Datenverarbeitung gebildet wird. Diese 
Aufhebung ist nicht exakt, da zwei Filter niemats genau 
gleich sind. Ebenso sind die Dopplerverschiebungen an 15 
unterschtedlichen Orten zu irgendeiner gegebenen Zeit 
unterschiedlich. Jedoch sind diese Resteffekte vernach- 
lassigbar, wie dies durch eine direkte Berechnung und 
durch aktuelle Experimente bestatigt wurde. 

In Fig. 5 ist ein Blockdiagramm der Digitalelektronik- 20 
einheit 37 dargestellt. Diese Digitalelektronikeinheit 37 
beinhaltet einen Signaiumsetzer 125» eine Korrelatoran- 
ordnung 127 mit einer Reihe von sieben identischen 
Korrelatoren, eine numerische Oszillatoranordnung 129 
mit einem entsprechenden Satz aus sieben identischen» 25 
numerischen Oszillatoren und eine Echtzeituhr 131, wo- 
bei die Korrelatoranordnung 127, die numerische OsziN 
latoranordnung 129 und die Echtzeituhr 131 iiber einen 
Datenbus 133 miteinander und mit dem Feldterminal- 
rechner 39 verbunden sind. Die erste Funktion des Si- 30 
gnalumsetzers 125 besteht darin, das analoge Obersei- 
tenbandsignal u(t), das analoge Unterseitenbandsignal 
l(t) und das analoge, sinusformige 5,115-MHz-Signal je- 
weits in ein binarwertiges "DigitaP- oder **Logik*-Signal 
umzusetzen, das zur Bearbeitung durch Qbliche Transi- 35 
stor-Logik-(TTL)-Schaltkreise geeignet ist 

Der Signaiumsetzer 125 erzeugt genau zwei Aus- 
gangsgrdfien. Eine dieser Ausgangsgrofien stellt eine 
binarwertige, TTL, periodische Rechteckwellenform 
mit einer Frequenz von 10,23 MHz dar, die durch Fre- 40 
quenzverdopplung der 5,115-MHz-Eingangsgrdfie er- 
zeugt wird. Dieses 10,23- MHz-Ausgangssignai dient als 
ein Taktsignal zur Zeitsteuerung aller digitalen Schal^ 
tungen. Dieses Taktsignal wird durch 1023 
( = 3x 11 x31) in der Echtzeituhr 131 dividiert, um einen 45 
Takt pro 100 Mikrosekunden zu erhalten. Weitere Divi- 
sionen mit aufeinanderfolgenden Faktoren von 10 er- 
zeugen dann eine vollstandige Dezimaldarsteliung der 
Zeit in Sekunden, wobei das niedrigstwertige Bit Einhei- 
ten von 10""* Sekunden reprasentiert. Die Zeit kann in 50 
dieser Form uber den Datenbus 133 immer abgelesen 
werden. Die Operationeh der Korrelatoranordniing 127, 
der numerischen Oszillatoranordnung 129 und des Feld- 
terminalrechners 39 werden alle durch die Echtzeituhr 
131 uber den Datenbus 133 geregelt. 55 

Die zweite "Digitar-Ausgangsgrdfle des Signalum- 
setzers 125 wird von den analogen u(t)- und I(t)-Ein- 
gangsgroBen abgeleitet und stellt eine binirwertige,. 
TTL-pegelige, unperiodische Wellenform dar. Diese 
AusgangsgroBe wird durch ein exklusives-TTL-NOR- eo 
Glied erzeugt, das zwei Eingange aufweist: ein Eingang 
verdeutlicht das Vorzeichen der u(t)-Eingangsgrd6e und 
der andere das Vorzeichen von l(t). Somit ist die Gatter- 
ausgangsgroBe "Wahr** (T oder binare 1) falls, und nur 
falls die analogen u(t)- und i(t)-Signale das gleiche Vbr- 65 
zeichen aufweisen. 

In Fig. 6 ist ein Blockdiagramm des Signalumsetzers 
125 dargestellt. Das analoge Singal u(t) stellt eine Ein- 
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gangsgroBe fur einen K^^^ator 135 dar. an dessen 
Ausgang ein TTL-Logikpegel erscheint, und zwar 
'Wahr", falls u(t) positiv ist und "Falsch" falls u(t) negativ 
ist. Dieses TTL-Logiksignal wird als eine EingangsgroBe 
fiir ein exkliisives TTL-NOR-Glied 137 verwendet. Das 
analoge Signal l(t) wird auf ahnliche Weise einem Kom- 
parator 139 zugefQhrt, dessen AusgangsgroBe als die 
andere EingangsgroBe des exklusiven NOR-GIieds 137 
verwendet wird Das sinusformige 5,1 15-MHz-Signal, 
das von dem Kristalloszillator 65 erhalten wird, stellt 
eine EingangsgroBe fur eine ubiiche analoge Frequenz- 
verdopplungsschaltung 141 dar, deren AusgangsgroBe 
einem dritten Komparator 143 zugefuhrt wird, um eine 
10,23*MHz-, rechteckwellenformige, TTL-pegelige Aus- 
gangsgrdBe zu erzeugen. Die 10,23-MHz- Ausgangsgro- 
Be wird ebenso als *Takt"- EingangsgroBe fiir ein Flip- 
Flop 145 verwendet. das die AusgangsgroBe des Gatters 
137 abtastet und halt; somit stellt die AusgangsgroBe 
des Flip-Flops 145 die exklusive NOR-Funktion der 
Vorzeichen von u(t) und l(t) dar. die mit einer gleichma- 
Bigen Geschwindigkeit von 10,23 x 10^ mal pro Sekunde 
abgetastet und wahrend der Abtastzeiten gehalten wird. 
Es ist im Stand der Technik der Funkinterferometrie 
allgemein bekannt, wie beschrieben z. B. von J. M. Mo- 
ran in einem Artikel in Methods of Experimental Phy- 
sics, voL 12, Teil C, Seiten 228—260, daB die binarwerti- 
ge Zeitfunktion U®L eine Fourier-Transformation oder 
-"Spektrum" aufweist, das eine gute Approximation so- 
wohl hinsichtlich der Phase als auch hinsichtlich der re- 
lativen Amplitude zu dem Fourier-Spektrum des analo- 
gen Produkts vori u(t)l(t) darstellt Die Genauigkeit der 
Approximation ist von den. Analogsignalen abhangig, 
die Zufallscharakter oder GauBschen Charakter aufwei- 
sen. Ebenso muB der Korrelationskoeffizient zwischen 
den beiden Eingangen in bezug auf seine GrdOe viel 
kleiner als 1 sein. (Praktisch "zittert" das Rauschen die 
Nichtlinearitaten der Komparatoren heraus. Das exklu- 
sive NOR-Glied 137 kann als eine Multiplizierschaltung 
angesehen werden, wobei jeder der Eingange Werte 
von -I- 1 und — 1 aufweist) Diese Bedingungen werden 
in dem vorliegenden System gut erfQllt Somit wird im 
folgenden der Logikpegel vom Flip-plop 145 einfach als 
Parstellung des Produkts u(t)l(t) betrachtet 

Das U^L-^Produkt** des Signalumsetzers 125 wird je- 
dem der sieben identischen Korrelatoren in der Korre- 
latoranordnung 1 27 parallel eingegeben. 

Ehe eine Beschriebung des Aufbaus der Korrelato- 
ranordnung 127 gegeben wird, werden seine Arbeits- 
prinzipien kurz erklart. 

In jedem Korrelator wird das u(t)l(t)-Produkt mit bi- 
naren Annaherungen an Sinus- und Cosinus-Zeitfunk- 
tionen korreliert, die durch einen entsprechenden der 
sieben numerischen Oszillatoren erzeugt werden. Die 
Frequenz des Oszillators wird durch deh Feldterminal- 
rechner 39 gemaB der Zeit, die von der Echtzeituhr 131 
angezeigt wird, gesteuert Zu irgendeinem gegebenen 
Zeitpunkt wird die Oszillatorfrequenz auf das zweifache 
der vorausgesagten Dopplerfrequenzverschiebung der 
1 575,42- MHz-Tragerwelle gesetzt, die von einem der 
Satelliten ubertragen wird. Jeweils ein Oszillator und ein 
Korrelator sind mit jedem der in Sicht befindlichen Sa- 
telliten; und zwar bis zu maximal sieben Satelliten, ver- 
bunden. (Falls mehr als sieben Satelliten immer in Sicht 
sind, konnen mehr numerische Oszillatoren und Korre- 
latoren in diesem System verwendet werden. In der Pra-. 
xis sind jedoch sieben ausreichend.) Falls die vorausge- 
sagte Dopplerverschiebung ausreichend nahe an die ak- 
tuelle Dopplerverschiebung heranreicht, w(erden dann 
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die AusgangsgroBen des Kori^j^s die Leistung und 
die Phase des Signals von dem speziellen Satelliten, fur 
den die Voraussage gemacht wurde, genau bestimmen 
und durch die Anwesenheit von Signalen von anderen 
SatelUten, die unterschiedliche Dopplerverschiebungen 
aufweisen, nicht nennenswert beeinfluBt 

Nachfolgend wird in mathematischen Ausdriicken die 
Arbeitsweise eines der numerischen Oszillatoren und 
seines zugehorigen Korrelators beschrieben: Als Funk- 
tion der Zeit t, die durch die Echtzeituhr 131 angezeigt 
wird, ist die vorausgesagte Dopplerfrequenzverschie- 
bung des Tragers des Satelliten durch den Ausdruck 
fp(t) gegeben. Der Wert von fp(t) wird aus der Tabelle 
aus vorberechneten Werten interpoliert, die fniher in 
dem Speicher des Feldterminalrechners gespeichert 
wurden, Der numerische Oszillator erzeugt zwei Zeit- 
funktionen: cos[2<Dp(t)] und sin[2 0p(t)], mit 90°-Pha- 
senverschiebung, wobei C>p(t) eine vorausgesagte Phase 
reprSsentiert, die eine Zeitfunktion darsiellt. Die Funk- 
tion Op(t) ist zu Anfang der Zeit to gleich Null, wenn der 
numerische Oszillator zu schwingen beginnt und an ir- 
gendeinem nachfolgenden Zeitpunkt Op(t) durch das In- 
tegral 



1 



*p(t) = 271 fp(t')dt' 



gegeben» wobet fp(t') den Augenblick von fp zu einer 
Zwischenzeit t' verdeutlicht Der Faktor von 2 n ist er- 
forderlich, da — wie ublich — die Frcquenz fp in Einhei- 
ten von Schwingungen pro Zeiteinheit und die Phase 4>p 
in Einheiten von Radian und hicht in Perioden gemessen 
wird. 

Nun bildet der [Correlator, der zwischen den Zeit- 
punkten to und ti arbeitet, GroBen a und b aus seinen 
Eingangswerten [u(t)l(t)l cos [2 <^p(t)] und sin [2 Op(t)], 
und zwar gemaB den Formeln 



u(l)l(t)cos(2<^p(t)ldt 



und 



• 1 

=1 



u(t)l(t)sinl2<l>p(t)]dt. 



iBe^ven 



uber das vorangegangen^B^sekundenintervall dar- 
stellt Diese AusgangsgroBe^verdeutlichen die Korre- 
iationen des Produkis u(t)l(i) mit den Kosinus- und Si- 
nusfunktionen. 

Wahrend des Einsekundenintervalls wird die Oszilla- 
torfrequenz fp(t) jede 0.1 Sekunde durch den Rechner 
aktualisiert, und zwar veranlaBt durch die 0,1-Sekun- 
dentakte der Echtzeituhr. Diese Aktualisierung ist er- 
forderlich^ dafi sich die Satelliten-Dopplerverschiebung 
andert, und zwar infolge der Bewegung der Satelliten 
relativ zum Feldterminal am Erdboden sowie der sich 
andemden Projektion der Relativgeschwindigkeit ent- 
lang der Sichtlinie mit einer Geschwindigkeit, die einen 
wesentlichen Bruchteil von 1 Hertz pro Sekunde dar- 
stellen kann. 

Nun konnen die KorrelatorausgangsgroBen a und b 
kombiniert werden. um Schatzungen der Leistung und 
der Tragerphase des Signals fur den speziellen Satelli- 
ten, fur den die Vorhersage fp(t) gemacht wurde, zu 
erhalten. 

Man definiert eine komplexe Zahl c. dessen Realteil 
gleich a urid dessen Imargin^rteil gleich b ist,d. h. 

c - a + jb, 

wobei j die Quadratwurzel von minus eins ist. Daraus 
folgt 

c a C < m* > < exp[2 j(0 — Op)] >, 

30 

wobei C ein positiver, reeller, konstanter Skalenfaktor. 
<m^> das Zeitmittel des Quadrats der GPS-Modula- 
tionsfunktion m(t) iiber das Integrationsintervall von to 
bis t| und <exp [2 j (<I>— <l>p)]> das Zeitmittel der kom- 

35 plexen Exponentialfunktion exp [2 j (<D— <l>p)] darstelU: 
Setzt man voraus, daB der Unterschied (O— Op) zwi- 
schen der empfangenen GPS-Tragersignalphase 
O ^ (l)(t) und der entsprechenden Voraussage Op = Op(t) 
sich nicht um einen wesentlichen Bruchteil einer 

40 Schwingung wahrend der Integrationszeit andert, dann 
ist die GroBe von c etwa proportional der durchschnitt- 
Hch empfangenen Leistung: 
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Das Integrationszeitintervall ti — to ist gleich einer Se- 
kunde, und die angezeigten Integrationen werden jede 
Sekunde durchgefuhrt Bei jedem Einsekundentakt sei- 
tens der Echtzeituhr werden die Werte der Integrale in 
die Speicherregister "strobiskopiert", die Intergrationen 
werden auf Null zurOckgesetzt, der numerische Oszilla- 
tor lauft wieder an. und eine neue Integrationsperiode 
beginnt. Somit Hefert der Korrelator am Ende jeder 
Zeitsekunde Gr6Ben a und b, die die Zeitmittel des Pro- 
dukts 

u(t)l(t)cos[2 0p(t)] 
bzw. des Produkts 
u(t)l(t)sin[2 0p{t)] 
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Id = (a2 + hY C < m2 > 

und der Winkel von c entspricht etwa zweimal der 
durchschnittlichen Phasendifferenz (O— Op): 

^c = tan"' (b/a) ^2 < (O-Op) >. 

Es ist zu bemerken. daB aus b und a der Winkel von c 
eindeutig modulo 2 n Radian bestimmt werden kann. 
Somit wird die Differenz (O— Op) modulo n Radian be- 
stimmt 

Damit die empfangene Signalleistung und die Trager- 
phase (modulo 7i) exakt aus a und b entspechend diesen 
Formeln bestimmt werden konnen, mussen zwei Bedin- 
gungen eingehalten werden: Erstens, wie bereits er- 
wahnt, muB sich die aktuelle Phase 0(t) von der voraus- 
gesagten Phase Op(t) um einen Betrag unterscheiden, 
der sich um viel weniger andert, als eine Schwingungs- 
periode wahrend der Integrationszeit von einer Sekun- 
de. Zweitens, der Rauschabsiand des Korrelatoraus- 
gangs, gegeben durch 

SNRc - (2/nXn/4XBeffTint)'^^F 

« (l/2XBeffTi„,)>^% 



- muB viel gro&er als Eins seir^^Bei Beff die effektive 
Bandbreite der Signale u(t) uncB^ist, etwa 5 x 10^ Hz; 
Tint stellt die Integrationszeit dar, gleich 1 Sekunde, und 
F stellt den Bruchteil der in u(t) und l(t) vorhandenen 
Leistung dar, der von dem GPS m(t)-Signal, und nicht 5 
vom Rauschen stammt Der Faktor (2/tc) erklart sich aus 
dem ICorrelatlonsverlust zwischen u(t) und l(t), der 
durch die Analog-Digital-Umwandlung dieser Signate 
durch die Komparatoren in dem Signalumsetzer her- 
vorgerufen wird. Der Faktor (71/4) erklart sich aus dem 10 
Veriust, der mit der Verwendung von' Rechteckwellen- 
approximationen der Sinus- und Kosinusfunktionen in 
dem Korrelator zusammenhangt. Die Wurzel aus dem 
Produkt von BeffTmt entspricht etwa 700. Demzufotge 
gilt die Relation: 15 

SNRc « 350 • F. 

Der Bruchteil F irgendeiner Seitenbandleistung, die 
von dem GPS-Satelliten abstammt. ist von dem Emp- 20 
fangsantennengewinn und der Rauschzahl des Emp- 
fangssystems abhangig. Fur die "MITES"- Antenna und 
das oben beschriebene Empfangssystem sowie fur einen 
Satellitenhdhenwinkel von etwa 20** ist aus einem Expe- 
riment bekannt, daB der Wert F etwa 0,03 Qbersteigt 25 
Demzufolge ergibt sich 

SNRc ^10, 

was fiir genaue Leistungs- und Phasenmessungen aus- 30 
reichend ist. Die Stand ardabweichung des Rauschens in 
jedem Teii der komptexen GroBe c, d. h. im Real- und 
Imagin^rteil, ist gegeben durch 
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ac ^ |c|/SNRc. 
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Die ersterwahnte Bedingung fQr die Genauigkeit der 
Messungen der Leistung und Phase, namlich daB 
(<D— (Dp) sich nicht um einen wesentlichen Bruchteil ei- 
ner Schwingung wahrend der Einsekundenintegrations- 40 
zeit andert. ist gleich der Bedingung, dafl der Unter- 
schied zwischen der aktuellen empfangenen Tragerfre- 
quenz f und der Oberlagerungsoszillatorfrequenz fo sich 
nicht von der vorausgesagten (numerischen Oszillator) 
Frequenz fp um einen wesentlichen Bruchteil von 1 45 
Hertz unterscheidet Dieser Bedingung wird in dem vor- 
liegenden System dadurch genugt, daB fiir die Frequenz 
des numerischen Oszillators eine selbstatige Regelung 
vorgesehen ist, um diese Frequenz nahe der aktuellen 
empfangenen Tragerfrequenz zu halten. Diese Rege- 50 
lung wird mit Hilfe eines einfachen Programms ausge- 
ubt, das durch den Feldterminalrechner 39 ausgefuhrt 
wird. Eine Beschreibung dieses Programms folgt 

Die komplexe Zahl c, die aus den a und b Korrelator- 
ausgSngeh am Ende des k-ten Einsekundenintegrations- 55 
intervalls gebildet wird, wird mit c(tk) gekennzeichnet, 
wobei tk die Zeit in der Mitte dieses Intervalls kenn- 
zeichnet. Zu der numerischen Oszillatorfrequenz fiir das 
(k -1-1)- Interval! wird ein Korrekturwert (corrective bi- 
as) 60 

K ♦ c [c(tk)c*(tk-t)]/27i Hertz 

hinzuaddiert, wobei K eine positive reelle FConstante 
kleiner als 1. < [ ] den Winkel der komplexen GroBe, die 65 
in den Klammern [ ] enthalten ist, und c*(tk- 1) die kom- 
plexe Konjugation der komplexen Zahl c des nachstvor- 
hergehenden (k — 1)-Intervalls darstellt. 



Das Arbeitsprtnzip did^^programms wird anhand 
des folgenden Beispiels ver^Rtlicht: Falls die Frequenz- 
vorhersage angenommen um 0,1 Hertz zu niedrig ist, 
dann wird der Winkel von c um 0,1 Schwingungsperiode 
in einer Sekunde voranschreiten, und die komplexe 
Grofle c(tk)c*(tk-i) wird einen Winkel von (+0,1) x(2n) 
Radian (plus geringes Rauschen mtt Mittelwert Null) 
aufweisen. Die Addition des Korrekturwerts, der in die- 
sem Falle positiv ist. wird die GroBe des negativen Feh- 
lers in der Frequenzvorhersage von (0,1 Hz) auf (1 — K) 
(0,1 Hz) verringern. 

Der Wert K muB groBer als 0 sein, oder keine Verrin- 
gerung eines Frequenzvorhersagefehlers wird von der 
Ruckkopplung resiiltieren. Der Wert muB kleiner als 1 
sein, Oder die ROckkopplung wird eine labile Schwin- 
gung bzw. Schwankung des Fehlers hervorrufen, und 
zwar infolge der Verzdgerung bei der Anwendung der 
FCorrektur. Der exakte Wert ist nicht kritisch, und der 
optimale Wert kanh durch Experiment bestimmt wer- 
den. In dem vorliegenden System wird ein Nominalwert 
von 0,5 verwendet. 

Ein wichtiger anderer Effekt dieser Frequenzriick- 
kopplung besteht darin, daB die numerische Oszillator- 
frequenz zu der aktuellen empfangjenen Tragerfrequenz 
'liingezogen" wird, und zwar von einer anf&nglichen 
Frequenz, die etwa mehrere Hertz oberhalb oder unter- 
halb liegen kann. Dieses "Hereinzieh^-Phanomen ist all- 
gemein bei der Technik der Phasen-oder Frequenznach- 
fiihrungsruckkopplungsschleifen allgemein bekannt, 
wie dies z. B. von Floyd M. Gardner in dem Buch "Pha- 
selock Techniques", veroffentlicht durch John Wiley & 
Sons, Inc^ New York, 1 966, diskutiert wird 

Die Bedeutung diieses "Hereinzieh^-Phinomens fur 
das vorliegende System, besteht darin. daB die a priori 
Kenntnis der Vermessungsmarkenposition keine gerin- 
gere Unsicherheit als einige Kilometer zu haben 
braucht 

Ein potentieller ungilnstiger Nebeneffekt dieses 
"Hereinzieh"-Ph2Lnomens besteht bei dem vorliegenden 
System darin, daB der numerische Oszillator der zur 
NachfQhrung eines speziellen Satelliten vorgesehen ist, 
statt dessen zu der Frequenz eines anderen Satelliten 
htngezogen wird, falls sich die Frequenz des letzteren in 
der Nahe des ersteren befmdet, und das Signal des letz- 
teren im Vergleich zum ersteren stark ist Um den Scha- 
den in Grenzen zu halten, der von solchen Vorfallen 
herriihren kann, enthalt das Feldterminal-Rechnerpro- 
gramm eine Vorschrift, die die GroBe des angewachse- 
nen Korrekturwerts, der zu der a priori Frequenzvor- 
hersage addiert werden kann, auf etwa 10 Hz begrenzt. 
Da der Unterschied zwischen zwei Satellitenfrequenzen 
sich andert, iiblicherweise etwa um 1 Hz pro Sekunde, 
folgt daraus. daB nur etwa MeBdaten von 10 Sekunden 
oder weniger als etwa 1 Prozent der gesamten, an einem 
MeBort empfangenen Daten durch Nachfuhrung eines 
falschen Satelliten hinfallig sein korinen. Die Erfahrung 
zeigt, daB dieser Prozentsatz unbeachtlich ist. 

in Fig. 7 ist ein Blockdiagramm eines Korrelatormo- 
duls 149 dargestellt, eines der sieben identischen Modu- 
le, die in der Korrelatoranordniing 127 vorhanden sind 
Alle sieben. Module haben die gleiche EingangsgroBe 
U®L, die dier U®L-AusgangsgrdBe des Signalumsetzers 
125 entspricht. Jedes Modul 149 empfangt ebenso eine 
"Kosinus"-Eingangsgr6Be und eine "Sinus"-Eingangs- 
groBe von einem entsprechenden Modul der sieben nu- 
merischen Oszillatormoduie. Die U^L- EingangsgroBe 
und die '*Kosinus*'-EingangsgroBe gelangen an eiti ex- 
klusives NOR-Glied 151, diessen AusgangsgrdBe die 
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EingangsgroQe eines **getakt^^HDigitaIzahlers 153 

darstellt. Die U®L-Eingangsgro^Knd die "Sinus^-Ein- , 
gangsgroBe gelangen an ein weiteres exklusives NOR- 
died 155, dessen AusgangsgroBe die Eingangsgr6Be 
eines weiteren Zahlers 157 darstellt. Einmal pro Sekun- 5 
de werden die Inhalte der Zahlerregister 153» 157 in 
entsprechenden Ausgangspuffern 159, 161 durch einen 
Impuls von der Echtzeituhr 131 in der Digitalelektronik- 
anorndung 37 verriegelt bzw. zwischengespeichert und 
die Zahler werden dann auf Null zuriickgesetzt Mit lo 
einem Takt von 10,23 MHz, gesteuert durch das 
*Takt"-Signal vom Signalumsetzer 125, wachst jeder 
Zahler 153, 157 urn eins, falls und nur falls seine Ein- 
gangsgroBe, die von seinem zugehdrigen exkiusiven 
NOR-Glied 151 brw. 155 abstammt, "Wahr" ist Somii 15 
zeigen am Ende jedes Einsekuhdenintervalls die Inhalte 
der Ausgangspuffer 159, 161 an, wie oft mal. und zwar 
zwischen Null und 10 230 OOOmal, die U®L- und die Ko- 
sinus/Sinus-Eingangsgroflen wahrend der vorgange- 
gangenen 1 Sekunde ubereinstimmten. Die Inhalte der 20 
Ausgangspuffer 159, 161 jedes Zahlers werden mit dem 
Datenbus 133 verbunden, uber den der Feldterminal- 
Rechner 39 zu jeder Sekunde die Inhalte liest. Jeder 
Zahler/Signalspeicher (Puffer) kann cine einzelne inte- 
grierte Schaltung darstellen, wie z. B. das 32-Bit- Bauele- 25 
ment, Modell Nr. LS7060, hergestellt von LSI Systems. 
Inc. 

Die GroBe a, die vorher durch die Kreuzmodulation 
zwischen [u(t)l(t)] und cos [2 <I>p(t)] defmiert wurde, wird 
in dem Feldterminal-Rechner 39 durch Subtraktion des 30 
Wertes 5 115 000 von der AusgangsgroBe des **Kosi- 
nus**-Zahlers und Division des Ergebnisses durch 
5 1 15 000 erhalten. Die GrbBe b wird ahnlich durch Sub- 
traktion des Wertes 5 115 000 von der AusgangsgrdBe 
des ''Sinus**- zahlers und Division des Ergebnisses durch 35 
5 1 15 000 erhalten. (Somit stellt die EinheitsgrdBe von a 
und b eine vollkommene Korrelation zwischen [u(t)l(t)] 
und der "Kosinus"- oder "Sinus"- Funktion dar. Ehe diese 
Ergebnisse im Speicher des Feldterminal-Rechners 39 
gespeichert werden, kann jede Zahl auf ebenso wenige 4 40 
Bit gekurzt warden, um Speicherplatz einzusparen.) 

In Fig. 8 ist eiii Blockdiagramm eines der sieben iden- 
tischen numerischen . Oszillatormodule 163 dargestellt, 
die in der numerischen Oszillatoranordnung 129 enthal- 
ten sind, wobei jedes numerische Oszillatormodul eine 45 
"Kosinus"- und eine "Sinus"-EingangsgrdBe zu einem 
Korrelatormodul 149 liefert. Jeder numerische Oszilla- 
tor 163 beinhaltet ein binares Phase«register 167 und 
ein binares Frequenzregister 169, einen binaren Addie- 
rer 171, ein exklusives NOR-Glied 173, einen Inverter 50 
175 und einen Frequenzteiler 177 

Das Phasenregister 167 und das Frequenzregister 169 
haben jeweils 32 Bit Der Addierer 171 ist ein 32-Bit-Ad- 
dierer. Die zu irgendeinem Zeitpunkt in dem Phasenre- 
gister 167 enthaltene Binarzahl stellt die Phase der Os- 55 
zillatorausgangsgroBe dar, wobei das hochstwertige Bit 
eine halbe Periode, das nachsthochstwertige Bit eine 
Vieriel Periode usw. darstellt. Die in dem Frequenzregi- 
ster 169 enthaltene Binarzahl stellt auf ahnliche Weise 
die Frequenz des Oszillators dar, wobei in diesem. Fall eo 
das hochstwertige Bit einen Wert von 155 000 Hz auf- 
weist, der gleich 1/66 der Periode des vom Signalumset- 
zer 125 stammenden lO,23-MHz-'Takt"-Signals ent- 
spricht. Der Addierer 171 summiert die in dem Fre- 
quenzregister 169 und dem Phasenregister 167 enthalte- 65 
nen Zahlen. Diese Summe wird in das Phasenregister 
167 geladen, unter Austausch der vorigen Inhalte, und 
zwar einmal pro Periode der AusgangsgroBe des Fre- 
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quenzteilers, der das tO,23-^^''Takt"- Signal durch ei- 
nen fesien Faktor von 33 t^^Das Phasenregister 167 
wird somit mit einer Geschwindigkeit von genau 
310 OOOmal pro Sekunde aktualisien. Der Betrag, um 
den die Phase nach jeder Aktualisierung fortschreitet, 
ist durch die inhalte des Frequenzregisters 169 gegeben. 
Das Frequenzregister 169 wird, wie bereits erwahnt, 
tOmal pro Sekunde uber den Datenbus 133 durch den 
Feldterminal-Rechner 39 aktualisiert. (Negative wie 
auch positive Frequenzen werden durch die Inhalte des 
Frequenzregisters dargestellt, und zwar unter Verwen- 
dung der bekannten Zweierkomplementmethode. Ge- 
maB dieser Obereinkunft wird das Negative einer Binar- 
zahl durch Komplementierung jedes Bits und anschlie- 
Bende Addition von eins gebildet Die groBte positive 
Zahl wird demgemslB dann dargestellt, falls das hochst- 
wertige Bit Null und alle anderen Bit eins sind. Falls das 
hochstwertige Bit eins ist, bedeutet dies, daB die Zahl 
negativ ist.) 

Die **Sinus"-AusgangsgroBe des numerischen Oszilla- 
tors 163 wird vom Inverter 175 erhalten, der das hochst- 
wertige Bit des Phasenregisters 167 invertiert. Die "Si- 
nus**- AusgangsgrdBe weist einen Wert von Eins auf, falls 
die Phase zwischen Null und einer positiven halben Pe- 
riode ist, und einen Wert von Null, falls die Phase zwi- 
schen einer halben und einer ganzen Periode ist (was 
das gleiche ist, wie wenn die Phase sich zwischen minus 
einer halben Periode und Null befindet). Die "ICosi- 
nus"- AusgangsgroBe des numerischen Oszillators 163 
wird von dem exkiusiven NOR-Glied 173 abgenommen, 
an dessen Eingangen das hochstwertige und das nachst- 
hochstwertige Bit des Phasenregisters 167 aniiegen. Die 
*'Kosinus''-Ausgangsgr6Be weist einen Wen von Eins 
auf, falls die Phase innerhalb von plus oder minus ein 
Viertel Periode, bezogen auf Null, ist. 

In Fig. 9 ist ein Blockdiagramm des Feldterminal- 
Rechners 39 gezeigt. Der Rechner 39 beinhaltet eine 
Zentraleinheit (CPU) 181. einen Programmspeicher 183, 
einen Datenspeicher 185, einen externen Zweirich- 
tungsdatenkanal 187, der mit einem Bedienungsterminal 
189 verbunden ist, und einen externen Zweirichtungsda- 
tenkanal 191, der mit einem Modulator-Demodulator 
(Modem) 193 verbunden ist, der weiterhin mit einer Te- 
lefonleitung, einem Funktelefon oder irgendeiner ande- 
ren Fernmeldeiibertragungseinrichtung 195 in Verbin- 
dung steht. Die Telle des Rechners 39 werden mit Hilfe 
eines Datenbusses untereinander verbunden, der auch 
zur Verbindung des Rechners 39 mit anderen Baugrup- 
pen des Feldterminals (vgl. Rg. 5) dient. 

Die Zentraleinheit CPU 181 kann z. B. ein Modell 
LSMl/2 (Teilenummer KDll-GC) der Digital Equip- 
ment Corporation (DEC) sein. Der Programmspeicher 
183 kann ein 32-lC-Byte programmierbarer Festwert- 
speicher sein. z. B. ein DEC, Teilenummer MRVll-C. 
Der Datenspeicher 185 kann ein 32-K-Byte Lese/ 
Schreibspeicher mit direktem Zugriff sein, z. B. ein DEC, 
Teilenummer MXVll-AC. Die beiden externen Zwei- 
richtungsdatenkanale (187 und 191) konnen die RS-232 
Seriendatenkanale sein, die in dem MXVl l-AC enthal- 
ten sind. Das Bedienungsterminal 189 kann das DEC 
Modell VT-100 oder irgendein gleichwertiges Serien- 
ASCIl-Terminal sein, das, wie das VT-100, mit der RS- 
232 Seriendatenschnittstelle des MXVll-AC oder 
durch irgendeine andere geeignete Datenkanaleinrich- 
tung mit dem Rechner verbunden werden kann. Der 
Modem 193 kann irgendeine, RS-232 kompatible. Stan- 
dardeinrichtung sein und kann vollstandig eliminieri 
werden, falls, wie erwahnt, der Feldterminal-Rechner 39 
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tions** in Methods of Ex perinj^M Physics, vol. 12, Teil 
C, Seiten 261 --276, 1976, besclWen,] 
Theoretisch ist die Grofie von S gegeben durch 

|Si| - C . G2(cos0j), 

wobei C eine Konstante und G den (Richtvermogen-) 
Gewinn einer Empfangsantenne darstellt, und zwar ge- 
. schrieben ais eine Funktion des Kosinus des Zenitwin- 
kels ©i des i-ten Satelliten. Man geht davon aus, da6 G 
unabhangig vom Azimut und derart normiert ist, daB die 
durch einen Kugelstrahler mit entsprechender Kreispo- 
larisation empfangene Leistung gleich 1 ist. Ftir den 
MITES-Antennenaufbau gilt 

G(cose) ==(U3) . (1 + cos 0)2 • sin2((3n/4)cos0) 
fQrO':S0^9O^ 
G(cpse)«0 . 
fur 90** ^8. 

Der Wert dieser Funktion ist etwa 2.46 am Zenit 
(0»O) und hat ein Maximum bei etwa 3,63 bei 0ss4O''» 
hat einen Einheitswert bei 0 a 72** und .erreicht 0, wenn 
090* erreicht. 

Der nachste Schritt bei der Bearbeitung der MeBda* 
ten, die von den beiden Interferometerterminals erhal- 
ten werden. besteht in der Summierung der komplexen 
Zahlen S|(t) iiber i; um eine Summe s(t) fur jeden MeB- 
zeitpunkt t zu erhalten: 



R(b) = i;iS(t,)| |S(t, 



S(t) = XSi(t). 



wobei die Summe sich iiber all die Satelliten erstreckt, 
die zu dem Zeitpunkt t beobachtet wurden. 

Der nachste Schritt bei der Bearbeitung der MeBda- 
ten besteht darin, einen Versuchswert 6 des Basisvek- 
tors F^aiiszuwahlen und von diesem Wert 6 eine Funk- 
tion der Zeit $(t) zu berechnen, die theoretisch den Wert 
repr^entiert, den S(t) gehabt hatten« falls der echte 
Wert b* des Basisvektors gleich dem Versuchswert 6 
ware: 



S(t) = XlAi(l)MBi(t)|exp(-j4iib-^(t)Ai|, 



wobei Xi die Funkweilenlange darstellt, die der empfan- 
genen Tragerfrequenz entspricht Das heiBt X\ = c/fj. Das 
Verfahren zur Auswahl eines Wertes von 6 wird unten 
beschrieben. Es ist zu benrierken, daB in dieser theoreti- 
schen Funktion §(t), im Gegensaiz zu der von Messun- 
gen abgeleiteten Funktion S(t), kein Term vorhanden ist, 
der die Oberlagerungsoszillatorphasendifferenz repra- 
sentiert Ebenso ist der konstante Skalenfaktor C weg- 
gelassen. 

AIs nachstes wird die GroBe von S(t) mit der Grofle 
von §(t) multipliziert, und das Produkt dieser GroBen 
wird uber alle MeBzeiten summiert, um einen Wert R(6) 
zu erhalten, der sowohl von 6 als auch naturlich von den 
Messungen abhangig ist: 



5 wobei ti die l-te der Reihe von etwa 5000 MeBzeiten 
darstellt. R(6) wird als "Mehrdeutigkeitsfunktion" be- 
zeichnet. 

Der nslchste Verarbeitungsschritt besteht darin, die 
Berechnung von R(f>) fur verschiedene Werte von 6 zu 
10 wiederholen und den speziellen Wert von 6 zu bestim- 
men. fur den die Funktion von R(B) den groBten Wert 
aufweist. Dieser Wert von B ist die gewunschte Bestim- 
mung des Basisvektors b! 

Der Versuchswert 6 des Basisvektors wird am An- 
15 fang^o gewahlt, daB er.der besten a priori Schatzung 
vpn b entspricht, die von einer unabhangigen Informa- 
tion der Lagen der Vermessungsmarken verfiigbar ist, 
wie z. B. die Lagen bzw. Orte, die man durch Identifizie- 
rung von Landmarken auf einer Karte erhalten kann. 
20 Die Maximierung von R(&) mit Bezug auf 6 wird durch 
Absuchen eines dreidimensionalen Raumes vorgenom- 
men, der an diesem Versuchswert 6 zentriert und groB 
genug ist, um die Unsicherheit der Anfangsschatzung 
einzuschlieBea Bei der Suche wird jeder Punkt eines 
25 gleichm^Big verteilten. dreidimensionalen Gitters ge- 
priift, um den einen Punkt zu lokalisieren, an dem R(B) 
maximal ist Der Gitterabstand betragt am Anfang 1 
Meter. Dann wird der Raum, der sich 2 Meter von die- 
sem einen Punkt des Maximums R(&) abgesucht. indem 
30 man ein Gitter mit einem Abstand von 20 Zentimetem 
priift. Das Maximum von k(6) wird in diesem feiner 
verteilten Gitter gefunden. Dann wird der Gitterab- 
stand und ebenso die lineare Ausdehnung des Gitters 
halbiert und die Suche wiederholt. Dieses Halbierungs- 
35 verfahren wird so lange durchgefuhrt, bis der Gitterab- 
stand unter 1 Millimeter liegt. Der Wert von 6, der 
schlieBlich R(6) maximiert, wird als gewunschte Bestim- 
mung des Basisvektors b herangezogen. Unter Verwen- 
dung einer Satellitenanzahl n gleich 5 kann eine Basis- 
40 vektorbestimmung mit Htlfe des erfindungsgemaBen 
Verfahrens mit einer Genauigkeit von etwa 5 Miilime- 
tern pro Koordinate fOr eine Basisl^nge von etwa 100 
. Metern vorgenommen werden. 

Das oben beschriebene Verfahren zur Verarbeitung 
45 der von zwei Interferometerterminals stammenden 
MeBdaten, um den Basisvektor zwischen den Terminals 
zu bestimmen. stellt eine Spezialisierung des allgemei- 
nen Verfahrens dar, das in der oben erwahnten Patent- 
anmeldung Nr. 3 05 142 beschrieben wird. Das in dieser 
50 Anmeldung offenbarte allgemeine Verfahren wird auch 
von Charles C Counselman und Sergei A. Gourevitch, 
in ''Miniature Interferometer Terminals for Earth Sur- 
veying: Ambiguity and Multipath with Global Positio- 
ning System**, veroffentlicht in IEEE Transactions on 
55 Geoscience an Remote Sensing, vol. GE-19., Nr, 4, Sei- 
ten 244—252, Oktober 1981 » beschrieben. 

Bei einem anderen Ausfuhrungsbeispiel des MeBda- 
tenverarbeitungsverfahrens gemaB der Erfindung wird 
auch eine Mehrdeutigkeitsfunktion R(6) aus den MeB- 
60 daiten und einem Versuchswert 6 der Basis gebildet; 
jedoch ist das Verfahren der Bildung der Funktion un- 
terschiedlich. In diesem Ausfuhrungsbeispiel, wie im 
vorhergehenden Ausfuhrungsbeispiel. wird die komple- 
xe Konjugierte von Ai(t) mit Bi(t) multipliziert, um ein. 
65 komplexes Produkt Si(t) zu erhalten : 



S.(t) = Ai^(t)Bi(t) 
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direkt mit dem Basisrechner IJ^^^unden ist. Der Da- 
tenbus 197 kann der LSI- 1 1 Q-^Kein, Die Echtzeituhr 
131» die numerische Oszillatoranordnung 129 und die 
Korrelatoranordnung 127 konnen mit dem Bus verbun- 
den werden, indem man sie auf Standardleiterplatten 5 
aufbaut, die direkt in den Plattenkantenverbinder der 
**Ruckebene" eines LSI-1 1 Rechnersystems eingesteckt 
werden. Derartige Leiterplatten stnd von DEC erhalt- 
lich. ausgestattet mit speziellen, integnerten Schaltun- 
gen« die die gesamte Datenkommunikation zwischen 10 
dem Q-Bus und den auf den Flatten aufgebauten spe- 
ziellen Interferometerterminal-Schaltkreisen durchfuh- 
ren konnen. 

Die in dem Speicher 185 des Feldterminal-Rechners 
39 gespeicherten Mefidaten weisen eine Zeitreihe aus 15 
komplexen Zahlen fur jeden der bis zu sieben beobach- 
teten Satelltten auf, wobei jede Sekunde eine derartige 
Zahl erhalten wird. Diese Daten werden fur eine Zeit- 
spanne von etwa 5000 Sekunden erworben, wahrendd- 
essen wenigstens zwei Satelliten standig beobachtet 20 
werden, wobei die Durchschnittszahl der beobachteten 
Satelliten wenigstens vier betragen soil. Fiir den i-ten 
Satelliten zur Zeit t wird die gegebene komplexe Gro&e 
mit Ai(t) bezeichnet, wobei die GrdQe dieser komplexen 
Zahl proportional der gemessenen Leistung des zu die- 25 
sem Zeitpunkt von diesen Satelliten empfangenen Si- 
gnals ist, wobei die Proportionalitatskonstante willkiir- 
lich, jedoch. fur alle Satelliten gleich ist. und wobei der 
Winkel der komplexen Zahl gleich dem Zweifachen der 
zu diesem Zeitpunkt fur den gleichen Satelliten gemes- 30 
senen Tragerphase ist, wobei die Phase jedes Satelliten 
auf das gleiche Oberlagerungs-Bezugsoszillatorsignal, 
namlich das 1 575,42- MHz-Signal bezogen wird, das 
durch die Oszillatoranordnung 57 des Feldterminals 
13-1 erzeugt wird. 35 

Die komplexen Daten Ai(t) i==l ..^ 7 werden durch 
den Feld terminal- Rechner 39 von den a und b Ausg^n- 
gen der sieben Korrelatoren 149 der Korrelatoranord- 
nung 127 wie fplgt abgeleitet: Fur den i-ten Korrelator 
gilt 40 



er^^H>i< 

w^ron 



Ai(t) »[a(t) + jb(t)]exp[2jOp(t)l 

wobei a(t) bzw. b(t) die normierte a bzw. b Ausgangs- 
groBe ftir das Einsekunden-'lntegrations"- oder Zahlin- 
tervall, zentriert beim Zeitpunkt t, darstellt; j ist die 
Wurzel aus minus Eins; und 2^p(i) ist das Zweifache der 
vorausgesagten Tragerphase des i-ten Satelliten zum 
Zeitpunkt t. Es ist zu bemerken. daQ die komplexe Zahl 
Ai(t) gleich der von dem i-ten Korrelatorausgang abge- 
leiteten komplexen Zahl c, multiplizien mit 
exp [2 j<l>p(t)] ist. Der Winkel von Ai stellt (zweimal) die 
empfangene Tragerphase bezogen auf (zweimal) die 
Phase des 1575,42-MHz-Uberlagerungsoszillatorsignals 
dar, wohingegen der Winkel von c auf (zweimal) die 
Summe der Bezugsoszillatorphase plus der Phase des 
numerischen Oszillators bezogen ist. 

Zum Zwecke der Erlauterung wird davon ausgegan- 
gen. da6 die Datenreihe |Ai(t)} von dem Feldterminal 
13-1 erzeugt wird, das sich am Anfang des Basisvektors 
befindet. Das andere Feldterminal 13-2, das das Feldter- 
minal am Ende des Basisvektors darstellt und die glei- 
chen Satelliten zu den gleichen Zeitpunkten wie das 
erste Terminal beobachtet, erzeugt Ai(t) entsprechende 
Daten. die mit Bi(t) gekennzeichnet werden. Die glei- 
chen Satelliten werden beobachtet, da beiden Terminals 
Vorhersagedaten vom gleichen Zentral- bzw. Basisrech- 
ner 15 zugefiihrt werden, der die Satelliten 1 bis 7 in 



genau einer Weise numer^^K>ie Beobachtungen an 
den beiden Terminals w^^n effektiv gleichzeitg 
durchgefuhrt, da die beiden Terminaluhren bzw. -takte 
unmittelbar vor den Beobachtungen synchronisiert 
wurden, und die Taktfrequenzen unterscheiden sich um 
einen unbedeutenden Betrag. (Der prinzipielle Effekt 
des Taktfrequenzunterschieds zwischen den Kristallos- 
zillatoren, die die Frequenzen der Takte liefern, besteht 
darin, daQ die Phasendifferenz zwischen den 
1575,42-MHz-Oszillatorbezugsfrequenzen sich andert) 
Es ist ohne Bedeutung. falls zu einer gegebenen Zeit ein 
spezieller Satellii von einem Terminal sichtbar, jedoch 
von einem anderen verborgen ist. In diesem Fall wird 
die GroBe von entweder Ai(t) oder Bi(t) einfach Null 
Oder in dessen Nahe sein. 

Die Operationen, die vom Zentralrechner 15 durch- 
gefOhrt werden. um die Bestimmung des Basisvektors 
des Interferometers abzuschlieBen, werden anschlie- 
Bend diskuliert. Hierzu werden dem Rechner die Lei- 
stungs- und PhasenmeBdaten zugefuhrt, die von den 
beiden an den Enden des Basisvektors angeordneten 
Feldterminals 13-1 und 13-2 gesammelt wurden. 

Der erste Schritt bei der Verarbeiiung der Ai(t)- und 
Bi(t)-Daten im Zentralrechner 15 besteht darin, die kom- 
plexe Konjugation von At(t), gekennzeichnet durch 
Ai*(t), mit Bi(t) zu multiplizieren. Das Produkt 

Si(t) = Ai»(t)Bi(t) 

hat einen Winkel ^Si(t), der gleich zweimal dem Unter- 
schied zwischen den gemessenen Phasen der Tragersi- 
gnale ist, die von den i-ten Satelliten an den beiden 
Terminals empfangen werden, wobei jede Phase mit 
Bezug auf den Ortsbezug- bzw. Oberlagerungsoszilla- 
tor in dem jeweiligen Terminal gemessen wurde. Dem- 
zufolge ist der Winkel von Si(t) init der Differenz zwi- 
schen den Phasen der Oberlagerungs- bzw. Ortsoszilla- 
toren und mit dem Basisvektor zwischen den Terminals 
durch die theoretische Relation 

^Si(t) ^ A^LO + (4nfi/c)r. sT(t) 



verknupft, wobei dieses A<Dlo die Phasendifferenz der 
OberlagerungsosziUatoren, fi die empfangene Frequenz 
fur den i-ten SateJIiten, etwa 1575,42 MHz^c die Licht- 

45 geschwindigkeit, b^der Basisvektor und Si(t) ein Ein- 
heitsvektor in Richlung des i-ten Satelliten, betrachtet 
zur Zeit t vom Mittelpunkt des Basisvektors aus, dar- 
stellt. (Diese Relation erzeugt den Winkel ^Si(t) in Ra- 
dian und nicht in Perioden. Da die Frequenz fi in Peri- 

50 oden bzw. Schwingungen pro Sekunde und nicht in Ra- 
dian bestimmt wird, muB ein Faktor 2n eingeschlossen 
werden. Der Grund dafiir, daB 4ji und nicht 2n hier 
erscheinen, ist darin zu sehen, daB jedes Feldterminal 
die empfangene Signalphase zweimal miBt.) Diese Rela- 

55 tion stellt eine Approximation insoweit dar, als diese 
eine Parallaxe zweiter Ordnung, Effekte des Obertra- 
gungsmediums, relativistische Mehrwegeffekte, Rau- 
schen usw., ignoriert Diese geringen Effekte werden 
hier aus Griinden der Klarheit vernachlassigt. Der mit 

60 der Vernachlassigung dieser Effekte einhergehende 
Fehler entspricht einem Basislinienfehler von weniger 
als etwa 1 cm bei einer Basislange von weniger als etwa 
1 km. [Mit Ausnahme des Rauscheffekts, der vdllig zu- 
fallig ist, ist es moglich, die oben vernachlassigien Effek- 

65 te nachzubilden, um eine genauere theoretische Darstel- 
lung von ^Si(t) zu erhalten. Diese Nachbildung wird z. B. 
in einem Artikel von I. I. Shapiro ''Estimation of astro- 
metric and geodetic parameters from VLBl observa- 



wobei Ai(t) eine komplexe Zaj^Bstellt, die reprasenta- 
tiv fur die Messungen des^^^^ils ist, das vom i-ten 
Sateiliten bei einem Interferometerterminal zur Zeit t 
empfangen wird, wobei die GroBe von Ai(t) proportio- 
nal der empfangenen Leistung ist, und der Winkel ^Ai(t) 
zweimal der Phase des Tragers relativ zum Oberlage- 
■ rungsoszillator des Terminals entspricht» und wobei Bi(t) 
gleich Ai(t) ist, mit der Ausnahme, daQ dieser Wert von 
dem anderen Terminal an dem anderen Ende des Basis- 
vektors abgeleitet wird. 

Als nachstes wird Si(t) mit einer gewissen komplexen 
Exponentialfunktion eines Versuchswerts 6 des Basis- 
vektors multipliziert, und das Produkt wird iiber alle zur 
Zeit t beobachteten Sateiliten summiert. urn eine Sum- 
me S(t) zu erhalten, die eine Funktion der Zeit und des 
Versuchswertes 6 darstellt: 
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und (d) nur ian^^Vin bezug auf die Phasen 
der TragerwellerWndern, die implizit in den 
empfangenen Signalen verhanden sind; und 
(e) zeitliches Mittein der Korrelationsproduk- 
te, wodurch Daten erhalten werden, die die 
Phase der implizit in dem von jedenri Sateiliten 
empfangenen Signal vorhandenen Tragerwel- 
le anzeigen» aus welbjien Phasen sich die Posi- 
tion des Erdpunktes berechnen liQt 

Hierzu 4 Seite(n) Zeichnungen 



S(t) = 2Si(t)expI-j4 7ibi;(t)/U 



20 



Wobei s i(t) einen Einheitsvekior in Richtung des i-ten 
Sateiliten zur Zeit t und X\ die Wellenlange des vom i-ten 
Sateiliten empfangenen Sisals darstellt (£s ist zu be- 
merken, daQ, falls b gleich b ist, der Winkel jedes Terms 
in der Summe uber i gleich A<I>lo ist, also unabhangig 
von i.) 

Als nachstes wird die GroBe von S(t) genommen und 
uber alle Beobachtungszeiten summiert, um die Funk- 
tion R(&) zu erhalten: 



25 



30 



R{b) - Si 5(1.) I, 



. * ■ . 35 

wobei ti die 1-te der etwa 5000 MeBzeiten darstellt. 

SchlieBlich wird der Wert von B, der R(6) maximiert, 
gefunden, und zwar durch den gleichen Suchvorgang, 
der in Verbindung mit dem urspriinglichen Datenverar- 
beitungsverfahren beschrieben wurde. Dieser Wert von^ 40 
6 stellt die gewunschte Bestimmung des Basisvektors b* 
dar. 

Dieses letztere Ausfuhrungsbetspiel ist rechnerisch 
effektiver als das zuerst beschriebene Ausfiihrungsbei- 
spiel. 45 

Patentanspruch 

Verfahren zum Bestimmen der Position eines Erd- 
punktes mit Hilfe modulierter Radiosignale, deren 50 
Trager unterdriickt ist, und die mit unbekannten 
pseudozufalligen Kodes moduliert sind- und von 
erdumkreisenden Sateiliten gesendet werden, ge- 
kennzeichnet durch folgende Schritte: 

(a) gleichzeitiges Empfangen der Signale von 55 
mehreren Sateiliten mit einer omnidirektiona- 
len Antenne am Erdpunkt, dessen Position zu 
biestimmen ist; 

(b) . Vorhersagen der Dopplerverschiebungen 
fiir die von der Antenne empfangenen Signale; eo 

(c) Korrelieren der empfangenen Signale mit 
den Vorhersagen fiir die Dopplerverschiebun- 
gen, um diese aufzuheben; 

(d) Korrelieren der von jeder Antenne emp- 
fangenen Signale mit sich selbst, um den unbe- 65 
kannten pseudozufalligen Modulationskode 

zu quadrieren, wodurch sich die Phasenwinkel 
der Korrelationsprodukte aus den Schritten (c) 
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